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Abstract. The concept of spin-orbital entanglement on superexchange bonds in tran- 
sition metal oxides is introduced and explained on several examples. It is shown that 
spin-orbital entanglement in superexchange models destabilizes the long-range (spin 
and orbital) order and may lead either to a disordered spin-liquid state or to novel 
phases at low temperature which arise from strongly frustrated interactions. Such novel 
ground states cannot be described within the conventionally used mean field theory 
which separates spin and orbital degrees of freedom. Even in cases where the ground 
states are disentangled, spin-orbital entanglement occurs in excited states and may 
become crucial for a correct description of physical properties at finite temperature. 
As an important example of this behaviour we present spin-orbital entanglement in the 
-RVO3 perovskites, with i?=La,Pr,- • -,Yb,Lu, where such finite temperature properties 
of these compounds can be understood only using entangled states: (i) thermal evo- 
lution of the optical spectral weights, (ii) the dependence of transition temperatures 
for the onset of orbital and magnetic order on the ionic radius in the phase diagram of 
the .RVO3 perovskites, and (Hi) dimerization observed in the magnon spectra for the 
C-type antiferromagnetic phase of YVO3. Finally, it is shown that joint spin-orbital 
excitations in an ordered phase with coexisting antiferromagnetic and alternating or- 
bital order introduces topological constraints for the hole propagation and will thus 
radically modify transport properties in doped Mott insulators where hole motion im- 
plies simultaneous spin and orbital excitations. 
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1. Introduction: Entanglement in many-body systems 

Superexchange models with spin-orbital entanglement on superexchange bonds, 
discovered by exact diagonalization of finite chains PQ, are a good recent example of 
entanglement in many-body systems. Entanglement is inherent to quantum mechanics 
and occurs in several systems. In general it means that quantum states have internal 
structure and cannot be represented as products of states which belong to different 
subspaces of the full Hilbert space [2HI]. This property of quantum states has gained 
renewed interest in recent years as it was found in several many-body quantum systems 
and it was realized that it may play a role in quantum information. It is shown below 
that it leads to measurable consequences in condensed matter systems with strongly 
correlated electrons, when orbital degrees of freedom are active. This new development 
concerns both model systems and the physical properties of Mott (or charge transfer) 
insulators — we summarize it shortly in the present topical review. 

Entanglement in quantum many-body systems is a broad field [3j and will not be 
discussed here as such. In the last few decades the interest in quantum entanglement 
has risen sharply in various formerly disconnected subfields of physics. At present 
these different communities come closer to each other in the search for universal ways 
of quantifying entanglement and developing algorithms to treat many-body quantum 
systems. An interested reader is encouraged to consult several review articles published 
recently on this subject — we name here only a few which focus on entanglement in: 
(i) many-body systems [5], (ii) quantum spin systems [6], (Hi) interacting fermionic 
and bosonic many-particle systems [7], {iv ) optical lattices [Sj, and finally, (t>) quantum 
cryptography and quantum communication [I]. Entanglement entropy plays a central 
role in these systems and is frequently used as a quantitative measure of entanglement [9]. 

Spin-orbital entanglement occurs either due to the relativistic on-site spin-orbit 
coupling or due to superexchange interactions on the bonds. While finite spin-orbit 
coupling introduces on-site entanglement, a qualitatively new and challenging situation 
is encountered when degenerate 3d orbitals of transition metal ions are partly filled and 
orbital degrees of freedom have to be treated on equal footing with electron spins in the 
effective spin-orbital superexchange model [10]. Here we mainly on the latter but some 
recent examples of entangled states in cases with strong spin-orbit coupling will also be 
mentioned for completeness at the end. 

When degenerate 3c? orbitals in a transition metal oxide are partly filled, realistic 
superexchange includes both orbital and spin degrees of freedom that are strongly 
interrelated [l0Tlll| . The microscopic models designed to describe realistic systems with 
strongly correlated and partly localized electrons include as well orbital interactions 
which follow from the orbital-lattice coupling and tune the orbital correlations. These 
latter interactions are rather strong in the e g orbital systems and stabilize the orbital 
order at rather high temperature, as for instance in LaMnOa [12J . It such cases it is well 
justified to treat the spin-orbital superexchange in the mean field (MF) approximation 
which separates orbital degrees of freedom and their dynamics from the spin ones [T3] . 
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The spin and orbital degrees of freedom order then in a complementary way and their 
order follows the classical Goodenough-Kanamori rules [14] . They were derived long ago 
from the microscopic insights concerning the structure of spin-orbital superexchange and 
predict that the antiferromagnetic (AF) order coexists with ferro-orbital (FO) order and 
ferromagnetic (FM) order coexists with alternating orbital (AO) order. Out of many 
examples which follow these rules, we mention here only the LaMnOs perovskite, with 
active e g orbitals and coexisting FM/AO order in ah planes and AF/FO order along the 
c axis [T2|[T5]. In this case spin and orbital operators indeed separate as the orbital order 
sets in at high temperature Too = 780 K and is already saturated when the A-type AF 
(A-AF) order occurs at Tn — 140 K. Therefore, the optical spectral weights measured 
in experiment are well described by the MF decoupling of spin and orbital degrees of 
freedom [16]. For this reason simple treatments of the models of manganites which use 
the MF approach [T71IT8] or Hartree-Fock decoupling [T9] are very successful in modeling 
the complex phase diagram of monolayer, bilayer and cubic manganites [20J. 

In a number of compounds with active orbital degrees of freedom where strong 
on-site Coulomb interactions localize electrons (or holes) and give rise to spin-orbital 
superexchange, two different types of long-range order compete with each other. A 
prominent example of this behaviour are the -RVO3 perovskites, where i?=Lu,Yb,- • -,La. 
In the case of perovskite vanadates the cubic symmetry is broken by GdFe03-like 
distortions and the xy orbitals are singly occupied at all V 3+ ions, while the second 
electron occupies the {yz, zx} doublet. One finds here two different AF phases for V 3+ 
ions in d 2 electronic configuration with S = 1 spin: (z) the C-type AF (C-AF) phase 
with staggered AF order in the ah planes accompanied by the FM order along the c axis, 
coexisting with weak AO order of active t 2g orbitals {yz,zx}, and (ii) the G-type AF 
(G-AF) phase with staggered AF order along all three cubic axes coexisting with robust 
C-type orbital order [21]. In fact, the coexisting orbital and magnetic order in these 
phases obey the Goodenough-Kanamori rules along the c axis. However, the situation 
in ab planes of this class of compounds is puzzling as two alternating orders coexist, 
both for spins and for orbitals. The reasons of this coexistence are more subtle — there 
also xy orbitals are singly occupied at every ion and they are in a FO state, while the 
spin order, driven mainly by them, is AF. So once again, the order in the ab planes can 
be understood using the Goodenough-Kanamori rules. However, the AF/AO order may 
be seen as entangled in the subspace of {yz, zx} orbitals and provides an interesting 
situation in doped systems, as we shall discuss below. 

Recent interest and theoretical progress in the understanding of spin-orbital 
superexchange models was triggered by the observation that orbital degeneracy 
significantly enhances quantum fluctuations which may even suppress long-range order 
when different types of symmetry broken states compete with each other near the 
quantum critical point [22]. The simplest and paradigmatic model in this context is 
the Kugel-Khomskii model introduced long ago for KCuF 3 [TO], a strongly correlated 
insulator with a single hole within degenerate e g orbitals at Cu 2+ ions in the d 9 electronic 
configuration. This model has two parameters which favour different types of symmetry 
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broken states: (i) Hund's exchange interaction, and (ii) the crystal field splitting of e g 
orbitals. When both these parameters are small the system is driven by its quantum 
nature — either long-range order disappears [22,23] or, for certain parameters, coexisting 
spin and orbital order might be stabilized by order out of disorder mechanism [21] . We 
shall discuss below to what extent the order is classical and show that spin-orbital 
entanglement manifests itself in the regime of most frustrated interactions. 

First, in this topical review we shall elucidate certain situations with the spin- 
orbital entanglement in the ground states (GSs). Such GSs are very challenging as 
there are no good methods in the theory to investigate them in a systematic way. It will 
become evident that quite different GSs arise, characterized by overestimated energy 
and incorrect correlation functions, when spin-orbital entanglement is neglected. 

Second, even when the GSs are not entangled, entanglement may be experimentally 
observed and has important consequences at finite temperature when the behaviour of 
the system is driven by low-energy excited states with spin-orbital entanglement. The 
existence of such states is a generic feature of any spin-orbital superexchange model and 
therefore the relevant question is only whether such states are accessible for thermal 
excitations. It will be shown that a rather exotic behaviour of the -RVO3 perovskites 
cannot be understood without including the spin-orbital entangled states. This point 
of view is supported by several experimental observations: (i) the thermal evolution of 
the optical spectral weights [25], (ii) the phase diagram of the -RVO3 perovskites [21], 
and (Hi) the observed dimerization in the magnon spectra of YVO3 [26J. 

An interesting situation arises also in doped Mott insulators, where doped holes 
introduce charge degrees of freedom which perturb the orbital order and frequently lead 
to phases with coexisting spin, charge and orbital order [27]. When orbital degrees of 
freedom are quenched, one finds that hole propagation occurs in the t-J model via a 
quasiparticle state that emerges due to quantum fluctuations in the spin background [28J . 
In t-J-like superexchange models with orbital degrees of freedom, hole propagation is 
either entirely suppressed by incoherent processes [29], or occurs by a rather subtle 
mechanism: either by off-diagonal orbital hopping in e g orbital systems [30J , or by next- 
nearest neighbour effective hopping in t 2g systems [3T1[32]. When both spin and orbital 
degrees of freedom may contribute, the situation is less clear as scattering on spin-flip 
processes introduces additional incoherence in hole propagation [33]. Surprisingly, it was 
realized only recently that spin-orbital entanglement introduces topological constraints 
for hole propagation in a Mott insulator with coexisting AF and AO order [Mj . and 
may thus have serious measurable consequences in doped i?V0 3 perovskites. 

The paper is organized as follows. In section 2 we explain general concepts of: 
(i) intrinsic frustration of orbital interactions, (ii) spin-orbital superexchange, and (Hi) 
its consequences for the magnetic exchange constants and the optical spectral weights. 
Next we describe spin-orbital entanglement in the GSs of spin-orbital models in section 
3. Such models are usually employed to explain the magnetic properties of transition 
metal oxides [13] and are also used to derive the optical spectral weights [35J . Therefore, 
even when the GSs are disentangled, entangled states have severe consequences on the 
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experimentally observed properties of some Mott insulators at finite temperature, and 
we describe in section 4 the properties of the perovskite -RVO3 systems as an example of 
such a complex behaviour driven by quantum entanglement. Spin-orbital entanglement 
may also occur in GSs in particular parameter regimes, and we provide two examples 
of this behaviour in section [5j (i) the Kugel-Khomskii model on a bilayer [36], and (ii) 
the d 1 model for t 2g electrons on a triangular lattice [371 138]. Entangled states may have 
also interesting consequences for hole propagation in ordered states — here we present 
the coupling of a hole to joint spin-orbital excitations [M] , see section 6. The paper is 
concluded in section 7, where a general discussion and main conclusions are presented. 

2. Orbital and spin-orbital superexchange 

2.1. Intrinsic frustration of orbital interactions 

Before presenting the spin-orbital entanglement, we first introduce the characteristic 
features of orbital superexchange interactions as obtained in case of spin polarized 
(FM) systems. These interactions are fundamentally different from spin superexchange 
which has high SU(2) symmetry and intrinsically frustrated. Frustration is one of the 
simplest concepts in physics with far reaching consequences [391 14"0]. The main and 
unusual feature of orbital interactions is their intrinsic frustration that follows from the 
directional nature of superexchange terms which contribute along different bonds and 
compete with one another [22]. This type of frustration does not follow from geometrical 
frustration and is best understood by considering a two-dimensional (2D) square lattice. 
In case of the Ising model frustration on a square lattice can be achieved, for instance, 
by changing signs of interactions along every second column and leaving the other 
interactions unchanged. In this case all plaquettes of the 2D lattice are frustrated as one 
of the interactions has the wrong sign but the spins order — the model is exactly solvable 
and has long-range order following the dominating interaction below a finite transition 
temperature [JT], being lower than the one of the 2D Ising model. We emphasize that 
this frustrated model is exactly solvable because it is still classical as the interactions 
concern only commuting {Sf} spin components. 

In contrast, the orbital interactions on a pseudocubic lattice are quantum, both for 
e g and t<i g orbitals, because they involve at least two pseudospin components [42]. Such 
models have different (typically cubic) symmetry from both the Z 2 symmetry of Ising 
and SU(2) symmetry of Heisenberg model, and are in general not exactly solvable on 
a 2D square lattice. We begin with the case of e g orbitals interacting within a 2D ab 
plane of K 2 CuF 4 compound; models for three-dimensional (3D) perovskites, for instance 
KCUF3 with d 9 electronic configurations of Cu 2+ ions or LaMn03 with d 4 configurations 
of Mn 3+ ions discussed below can be easily obtained as a straightforward generalization 
of the 2D model using the cubic symmetry of orbital interactions [43]. Two e g orbital 
states, 
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are the eigenstates of the r[ = \{n iz — n ix ) orbital operator for pseudospin r = |, 
where {ni Z ,n ix } are hole number operators at site i. 

The origin of intrinsic frustration in the e g orbital superexchange is best realized 
by considering its form [43 J, 

H eg = J orb J2 U l] ^ ] - \) , (2) 



T (o.6) = 

1 " 4 



where the bond is oriented along one of the cubic axes 7 = a, b, c [43] . Here the orbital 
interaction J or b > follows from the energy of the high-spin charge excitation [13] . The 
pseudospin operators take a different form depending on the bond direction and are 
defined as follows, 

1 ' -of ± V3of) , rt ] = \al (3) 

where a* are Pauli matrices and the sign ± in is selected for a bond (ij) along 
a (b) axis. Thus for the ab plane one has two linear combinations of {erf, erf} Pauli 
matrices, and these interactions favour AOs on each bond, being the eigenstates of the o~f 
Pauli matrix as the interactions oc cf <rj are here the strongest ones. We emphasize that 
the interactions in equation (j2j) are fundamentally different from the SU(2)-symmetric 
spin interactions, as they: (i) obey only lower cubic symmetry, (ii) are Ising-like and 
select only one component of the pseudospin interaction along each bond which favours 
pairs of orthogonal orbitals, i.e., oriented along the bond (z-like) and the orthogonal 
to it lying in the plane perpendicular to the bond (x-like). This manifests the intrinsic 
frustration of orbital interactions in the e g orbital case [22J. In fact, the interactions 



in equation (JSJ) are Ising-like and classical only in the one- dimensional (ID) model 
but in general they are not. However, due to the gap which opens in orbital excitations 
in the 2D model, the quantum corrections generated by them are rather small [43J. 

The essence of orbital frustration which characterizes the e g orbital superexchange 
fl2]) is captured by the 2D compass model [45] which arises by increasing frustration 
from the 2D e g orbital model to the maximal frustration [46]. One considers then the 
2D model (or an exactly solvable model compass ladder [47] ) that interpolates between 
the classical Ising one and the compass one passing through the e g orbital model. The 
orbital interactions J or b are equal along both rows and columns but select two orthogonal 
pseudospin components [451148] : 

h 2d = jJ2 t ^ + j J2 t ^- ( 4 ) 

Usually one considers the AF case (J > 0) but the FM model (J < 0) is equivalent 
and equally frustrated. Intersite interactions in the 2D compass model are descibed by 
products T^Tj 1 of pseudospin components with a — x,y, z, 

T i = ' T i = 2 l ' r< = 2 ' ' 

rather than by pseudospin scalar products T{ • Tj. As explained below, such scalar 
products arise for the superexchange interactions with active t 2g orbitals degrees of 
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freedom which allow hopping processes for a pair of them in each 2D plane in the cubic 
system. Instead in the compass model (jlj) the ifr? interactions for bonds (ij) along the 
a axis compete with the rfr? ones along the b axis [48J. Also the ID compass model 
with alternating rf 2zr| i+1 and rf^r^ interactions [H] is intrinsically frustrated. 

Recently the 2D compass model was investigated by Monte Carlo simulations and 
the existence of a phase transition at finite temperature was established [50]. The 
ordered GS is degenerate and its different states correspond to either eigenstates of 
{rf} pseudospin components ordered along the rows, or eigenstates of pseudospin 
{rf} components ordered along the columns [51]. Although this GS is destabilized 
by infinitesimal pseudospin interactions, the structure of the lowest excited states stays 
unchanged and corresponds to the flips of spin columns [52] . It has been suggested that 
the model could serve as an effective model for protected qubits and such states realized 
by Josephson arrays [53] could play a role in quantum communication. Indeed, first 
experimental successes in constructing special networks of Josephson junctions that are 
designed following the compass model were reported recently [5%] . 

Although frustration still increases by going from the 2D to 3D e g orbital model, 
there are indications from recent Monte Carlo simulations that the GS is ordered [55,56j. 
Disorder occurs here by doping which leads to the 3D orbital liquid state [57] that plays 
a prominent role in the FM metallic manganites and provides an explanation for the 
observed magnon dispersion [HE]- The case of the 3D compass model is more subtle. It 
was concluded from the high temperature expansion that the ordered state is excluded 
here [59], but this result was challenged recently [55] and further studies are needed 
to establish whether this model could indeed serve as an example of an orbital liquid 
phase. 

The orbital models for t<i g orbitals contain more terms and were less studied up to 
now. The generic form contains scalar products of two pseudospins r = 1/2 along each 
direction in the cubic lattice, defined as in equations ([5]) [60J, 




where the operators r/ = {Tf,rf,rf}^ depend on the bond direction, i.e., the 
pseudospin components are defined here in different subspaces of the Hilbert space, 
depending on the pair of active t2 9 orbitals. This form follows from the fact that only 
two t2 g orbitals allow for electron hopping along a given cubic direction 7, while the third 
orbital is inactive, see also below and section 4. For finite Hund's exchange additional 
terms arise and the orbital state is disordered [6Tj . A stable orbital ordered state was 
found here in the 3D orbital model for YTi03, when the spins decouple from ti g orbitals 
in the FM phase |62j . 

2.2. Spin-orbital superexchange models 

In transition metal compounds with large on-site Coulomb interactions charge 
fluctuations are suppressed and electrons partly localize. This happens when intraorbital 
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Coulomb interaction U is large compared to the effective d — d hopping element t, 
where t is either the (ddcr) or (ddn) effective d — d hopping element for e g and t 2g 
systems, respectively, that arises via hybridization with ligand orbitals. In the regime 
of t <C U, correlated (Mott or charge-transfer) insulators arise in undoped compounds, 
while doping leads to interesting phenomena in strongly correlated electron systems 
with charge fluctuations only between two neighbouring electronic configurations [63J. 
Intersite charge excitations may be then treated within perturbation theory, while the 
hopping processes which do not cost high local Coulomb energy U and occur in doped 
systems are treated in the restricted Hilbert space. A well known example of this 
description is the t-J model, used widely to describe the physical properties of high-T c 
superconductors, but derived one decade before their discovery [64]. 

Here we concentrate on Mott insulators with transition metal ions in d n electronic 
configuration and active orbital degrees of freedom, where the effective low-energy 
Hamiltonians contain spin-orbital superexchange, described within spin-orbital models 
pHHH]. Such models are derived using the realistic multiplet structure of the excited 
states of transition metal ions which arise in d™ l d r J l =± <i- m+1 ^<i ( ; m ~ 1 ^ intersite charge 
excitations. As the multiplet structure depends on the electron number m [65], with 
some examples given in [13], these models are specific to a given family of compounds. 
As a representative example we consider here in more detail the case of the -RVO3 
perovskites (see section 4), with 5=1 spin stabilized by Hund's exchange and active 
t 2g orbitals at V 3+ ions in d 2 (m = 2) electronic configuration. In a cubic perovskite all 
three t 2g orbitals are degenerate and the kinetic energy electrons is given by: 

H t = ~ t Yl Yl ( d laJ jaa + d) au d ia ^j , (7) 

7 01(7)7° 

where d\ aa is electron creation operator for an electron at site i in orbital state a with 
spin o =t, 4- The summation runs over three cubic axes, 7 = a, b, c, and involves the 
bonds (ij) || 7 along them. The hopping t results from transitions via an intermediate 
0(2p n ) orbital and conserves the active t 2g orbital flavor a (7). Thus the hopping t is an 
effective (ddir) element that originates from two subsequent d—p hopping processes along 
each V— O— V bond in the -RVO3 perovskite structure. For e g systems the derivation is 
similar but the effective hopping elements follow from the hybridization with 0(2p a ) 
orbitals. The effective hopping follows from the charge-transfer model with p — d 
hybridization t p ^ and charge-transfer energy A [66J, and one expects in the present 
vanadate case t = t"i d /A ~ 0.2 eV [57] . Only two t 2g orbitals, labelled by 0(7), are active 
along each bond (y)||7 and contribute to the kinetic energy ([7]), while the third one lies 
in the plane perpendicular to the 7 axis and the hopping via the intermediate oxygen 2p n 
(or 2p u ) oxygen orbital is forbidden by symmetry. This justifies a simplified notation 
used below, with the orbitals defined by the axis direction which is perpendicular to 
their plane: 

\ a ) = \yz), \ b ) = \ zx ), l c ) = M- ( 8 ) 

In this case only t 2g orbitals are partly filled by electrons, and it suffices to consider 
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local Coulomb interactions between t 2g electrons at V 3+ ions described by the degenerate 
Hubbard Hamiltonian [68J, 



Hint = U ^2n iat n ial + ( U - ^Jh) 

ia ^ ' i,a<[ 



i,a</3 



i,a<0 

— 2Jh Si a ■ Si/3 . (9) 

i,a</3 

Here Ui a = Ui aa is the electron density operator in orbital a = a, b, c at site i, and 
spin operators Si a = {Sf a , Sf a , S? a } for orbital a at site i are related to fermion operators 
in the standard way, i.e., 

^ia = d io ^d ia ± , Sf a = -{n ia ^ — n ia i) . (10) 

The first term in ([9]) describes the largest intraorbital Coulomb interaction U for a pair of 
electrons with antiparallel spins in orbital a. The second term stands for the interorbital 
Coulomb (density) interaction, the third one is called frequently the "pair-hopping" 
term, and the last one is Hund's exchange Jh- The choice of coefficients in (jUJ) guarantees 
that the interactions satisfy the rotational invariance in the orbital space [HE]. This 
Hamiltonian is exact when it describes only one representation of the cubic symmetry 
group (here t 2ff orbitals which are partly occupied in the cubic vanadates) — then the 
on-site interactions are given by two parameters: (i) the intraorbital Coulomb element 
U, and (ii) Hund's exchange element Jh- These elements may be expressed by the 
Racah parameters {A, B, C}. For t 25 electrons considered in section 4 one finds [131165] : 

U =A + 4B + 3C, (11) 
J H = W + C. (12) 

Hund's exchange (and interorbital Coulomb interaction) is in general anisotropic and 
depends on the pair of interacting orbital states. For instance, the corresponding e g 
Hund's exchange element is Jh = + C. More details are given in |13j . 

In the limit of large U (t U), the effective low-energy spin-orbital superexchange 
interactions arise by considering all the contributions which originate from possible 
virtual excitations d™d™ ^ d™ +1 d^ 1 " 1 . The general structure of spin-orbital 
superexchange [13], 

Hj = JJ2{^ (3 • 3 + S 2 ) + ^} , (13) 

—* —* 

involves interactions between SU(2)-symmetric spin scalar products Si ■ Sj on each bond 
(ij), connecting two nearest-neighbor transition metal ions, each one coupled to orbital 
operators {-fj, fj}. The orbital operators are given in section [27T1 and obey only much 
lower symmetry (at most cubic for a cubic lattice). These operators contribute to the 
form of orbital operators and which depend on the model. They involve the 
active orbitals on each bond (ij) along direction 7. 
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In order to derive magnetic excitations for the systems with orbital degeneracy 
one usually derives magnetic exchange constants for a bond (ij) by averaging over the 
orbital operators in equation ( {TBI using the GS |$o) with decoupled spin and orbital 
operators, 

^ = <*o|^ 7) |<&o>. (14) 

This procedure assumes implicitly that spin and orbital operators can be decoupled 
from each other in the MF approach and ignores the possibility of entanglement and 
composite spin-orbital excitations introduced in [22]. Inter alia, such excitations play a 
prominent role in destabilizing the classical AF long-range order in the d 9 spin-orbital 
model [23] . 

The energy scale for the superexchange is given by 

J=l- (15) 

where t is the relevant effective hopping element and U is the intraorbital Coulomb 
element defined in Q. As several charge excitations contribute to the superexchange 
(fl3|) . the balance between competing terms depends on Hund's exchange, namely on 

V=% (16) 

This is the only parameter which decides about the strength of particular interactions in 
the superexchange and finally also about the type of magnetic correlations or symmetry 
breaking in the GS obtained at orbital degeneracy and favoured by the model. 

We would like to emphasize here that the same charge excitations which decide 
about the spin-orbital superexchange ffl3|) contribute as well to the optical conductivity. 
In this case they appear at distinct energies of individual charge excitations that occur at 
a given ionic filling, and depend on the multiplet structure of the excited states arising 
from intersite charge transitions [35J. Each of these excitations involves a different 
intermediate state in the multiplet structure of at least one of the ions after the charge 
excitation, i.e., either in the d m+1 or in the d m ~ l configuration or in both, depending 
on the actual process and on the value of m [35]. This feature made it possible to 
relate the averages of these different excitations to the spectral weights in the optical 
spectroscopy [35], and this principle serves now as a theoretical tool to analyse and 
explain the observed anisotropy and temperature dependence of the spectral weights in 
the optical spectra [13J. 

As the charge excitations correspond to particular expressions in the spin-orbital 
space, it is important to rewrite ther superexchange Hamiltonian ( fl3|l by decomposing it 
into individual terms on each bond (ij) that stem from particular excited states labeled 
by n [35], 

Here the superexchange constant (fT5|) was included in the operators Hn (ij). As 
explained above, the spectral weight in the optical spectroscopy is given in a 
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correlated insulator by the same virtual charge excitations that contribute also to the 
superexchange. They define the individual kinetic energy terms Kn^ along the axis 
7, which can be determined from the superexchange ( {TBI using the Hellman-Feynman 
theorem [69] , 



For convenience, we define them here as positive quantities, Kn > 0. 

The spectral weights (ITS1) are defined by the bond correlation functions and their 
changes with increasing temperature which decide about the temperature dependence of 
the optical spectrum. To describe experimental observations it is therefore important to 
analyse the various multiplet excitations separately, as they depend on these correlations 
in a different way, and will also contribute to a quite different temperature dependence, 
as we show in this paper on the example of LaVC>3 in section 4. Such an analysis is 
of course possible in each the respective spin-orbital superexchange models are 

derived by considering all different types of charge excitations individually [T3] . 

In some cases, however, spin dynamics separates from the orbital one in the GS and 
MF factorization of spin and orbital operators is indeed allowed. This happens when the 
orbital order is stabilized to a large extent by strong interactions with the lattice which 
undergoes Jahn- Teller distortions. A good example of this behaviour is LaMnOs, where 
the superexchange and the Jahn- Teller effect help each other [12] in stabilizing the orbital 
order which occurs below the relatively high transition temperature Too — 780 K. The 
A-AF spin order is observed below a much lower Neel temperature Tn — 140 K [16J. In 
this case a very satisfactory description of the experimental results for the spectral 
weight distribution in the optical spectroscopy is obtained using disentangled spin- 
orbital superexchange, both for the high spin [13] and low spin [16] optical excitations. 
Below we focus on some examples of a more complex behaviour driven by the spin-orbital 
entanglement. 

3. Spin-orbital entanglement 

Before presenting the essence of spin-orbital entanglement, we would like to remind the 
reader severe consequences of the MF approximation used frequently to investigate spin 
and/or orbital models. In this approach quantum fluctuations are neglected and only 
qualitative conclusions concerning possible symmetry breaking can be drawn. In spin- 
orbital systems the MF approach may be applied in different ways, and here we focus 
on the decoupling of spin and orbital degrees of freedom presented below for a single 
bond in section 13.11 When only spin-orbital coupling is treated in the MF approach, 
spin and orbital operators are disentangled. We show in section [3721 and section [3731 that 
this approximation is unsatisfactory in many situations. 

To detect spin-orbital entanglement in the GS we evaluate intersite spin, orbital and 
joint spin-orbital bond correlations in several models, defined as follows for a nearest- 
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neighbour bond (ij) (we keep the notation general here, for the ID chain j = i + 1) [Tj: 



The above general expressions imply averaging over the exact GS found from Lanczos 



the models for t 2g orbitals investigated in section 13.31 While Sij and T y - correlations 
indicate the tendency towards particular spin and orbital order, C y - quantifies the spin- 
orbital entanglement — if C y - < spin and orbital operators are entangled and the 
MF approximation, i.e., decoupling of spin and pseudospin operators in ( Fl3l) . cannot be 
applied as it generates uncontrollable errors. 

3.1. Exact versus mean field states for a bond 

The spin-orbital superexchange ( 1T3"]) takes the simplest form when Hund's exchange is 
absent (Jjj = 0) — then the orbital operators J-^ which couple to the scalar spin 
product (Si ■ Sj) are particularly simple. In case of e g orbitals they are just projection 
operators on the active directional orbital [70], while for t 2g orbitals they give a scalar 
product of r = ~ pseudospins which represent two active orbital flavours along the 
bond direction (171) . As an example, we consider first a bond in the ID SU(2)®SU(2) 
model for S — 1 spins of V 3+ ions coupled to r = | pseudospins. Here we show that 
quantum fluctuations in both spin and orbital system are crucial to reproduce faithfully 
the energy spectrum of a bond with the spin-orbital superexchange, 



Such superexchange interactions are obtained in a perovskite RV0 5 vanadate along the 
c axis in absence of Hund's exchange (rj = 0) [71]. The above form is convenient for 
further discussion and we neglected here a constant term (which would play a role for 
the optical spectral weights [35J). 

While the energy spectrum for a single bond can be easily solved exactly using 
the SU(2)®SU(2) symmetry and the classification of quantum states by the total spin 
Sf = Si + 5*2 and orbital (pseudospin) r t = fi + f 2 operators, in a solid or in an infinite 
chain this is not the case and one has to employ some approximation. Frequently the MF 
theory is used which has severe limitations as the quantum fluctuations on the bonds 
are then neglected. This can be seen by solving the problem of a single bond ( 1221 . 
either by replacing the scalar products of both spin and orbital operators by their zth 
components (T-L^^i or by treating in the MF approximation only the (more classical) 
spin scalar product for S = 1 spins (Hf 2 ): 




(19) 
(20) 
(21) 



diagonalization of a finite cluster and are valid for S 



| and S = 1 encountered in 




(22) 




(23) 
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Figure 1. Eigenenergies and their degeneracy (numbers above vertical lines) for 
a single spin-orbital superexchange bond in RVO3 at 77 = 0, described by: (a) 
classical Ising-like interactions ([23| . (b) Ising-like interactions for S = 1 spins but 
full quantum interaction for r = | pseudospins and (c) the full quantum model 
(|22|) . Spin-orbital entanglement is absent here, but orbital fluctuations change the 
energy spectrum in a dramatic way and stabilize the high-spin state. 

W 12 =\J(SISZ + 1) (n.r 2 + ^). (24) 

Note that the Ising-like Hamiltonian ( 123]) is nonnegative by construction, so the GS 
energy E = is obtained when at least one of the following conditions is satisfied: 
either SfS^ = —1 or rfr| = — |. This property gives a rather high degeneracy d = 22 
of the GS, see figure [H(a). At the same time the highest energy in the spectrum with 
degeneracy d — 4 is accurately reproduced when SfS^ = +1 and rfrf = +7, i.e., for four 
possible configurations where quantum fluctuations do not contribute (2 configurations 
for spin and 2 for orbital operators). Note that although the highest energy E — | J is 
correct, the degeneracy of this state at is not correctly reproduced. 

The full bond Hamiltonian (|22l) has however also negative eigenvalues and the 
energy spectrum which includes quantum fluctuations starts at a much lower energy 
Ei = —\J, see figure dc), obtained for a pseudospin singlet r t = and the high spin 
S t = 2 state. The exact spectrum, shown in figure [U^c), is obtained by considering all 
eigenstates with S — 0, 1, 2 and r D = 0, 1. In fact, it is even sufficient here to consider 
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only the exact values of the f\ ■ r 2 operator for the orbital singlet (r D = 0) and orbital 
triplet (t = 1) states and to keep the Ising form of the spin interaction (123]) to reproduce 
all exact eigenenergies for a single bond, see figure E(b). However, only when the spin 
configurations and the eigenstates of total spin S are used, the degeneracies of all the 
states in the spectrum are incorrect, cf. figures [T^b) and[H(c). Therefore, we note that 
the thermodynamic properties determined in the MF theory are not free from systematic 
errors. We show below that full spin and orbital dynamics on the bond plays a very 
important role in systems where orbital correlations contribute by quantum effects. 

Yet, there is even a more serious problem which concerns the spin-orbital 
superexchange ( fl~3l) the possible entanglement of spin and orbital operators. 
Entanglement means here that wave functions cannot be written as products of spin 
and orbital states, similar to entangled spin singlet wave function |2H1], where it is not 
just one eigenstate of the S- spin operator at site i. As the wave functions in the present 
example obey the SU(2) symmetry in spin and orbital subspace, entanglement does not 
occur for a single bond, where only individual spin singlet or orbital singlet states are 
entangled. However, spin-orbital entanglement is a characteristic feature of any spin- 
orbital model in a larger system, both on a finite cluster and in the thermodynamic 
limit. There the same spin and orbital operators participate in interactions along several 
bonds, and we show below that this is also the case for the interaction in equation ([22]) . 
The essence of this type of entanglement is explained in the following section. 

3.2. Entanglement in the SU(2)®SU(2) spin-orbital model 

Before demonstrating spin-orbital entanglement in more realistic models which apply to 
systems with active t 2g orbital degrees of freedom, we consider first the SU(2)®SU(2) 
spin-orbital model, with S = ~ spins and orbital interactions given by a scalar product 
Ti ■ Tj for r = | pseudospin operators. Two components of r = | pseudospin stand for 
two active t 2g orbitals along the 7 axis in a cubic (perovskite) lattice. Indeed, for a 
given cubic direction only two t 2g orbitals are active [66J and the orbital superexchange 
has SU(2) symmetry when Hund's exchange is neglected. 

Here we shall analyse the energy in the ID SU(2)xSU(2) model, 

n j = J £ (3 ■ &n + x ) { 7i ■ + y) ' ( 25 ) 

i 

for spins S = \ and pseudospins r = ~, using exact diagonalization of finite chains 
with periodic boundary conditions (PBC). The model (125]) has two parameters x and y. 
Its phase diagram in the (x, y) plane consists of five distinct phases which result from 
the competition between effective AF and FM spin, as well as AO and FO pseudospin 
exchange interactions on the bonds [72] . First of all, the spin and pseudospin correlations 
are FM/FO, if x < —\ and y < — |. Then the GS is characterized by the maximal 
values of both total quantum numbers, S t = % = N/2, where N is the chain length; its 
degeneracy is d = (N + l) 2 . Two other phases have either FM spin or FO pseudospin 
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Figure 2. Energy spectrum for the spin-orbital chain ([25]) of N = 4 sites with 
periodic boundaru conditions and x = y = p, for increasing p. At p = — \ the GS 
changes from high spin-orbital state (S = T = 2 with degeneracy d = 25) to spin- 
orbital singlet state (S — T = with d = 1). The point with SU(4) symmetry (p = t 
with g? = 1) and the special SOVB point (p = | with <i = 2) are marked by vertical 
dashed lines. This figure is reproduced from |76j . 

order, accompanied by alternating order in the other channel, i.e., AF order for the FO 
phase and AO order for the FM phase. 

A unique feature of FM state is that it is an eigenstate of the Heisenberg exchange 
operator. The applies to the orbital interactions, so the quantum fluctuations are 
entirely suppressed in the FM/FO phase. They are also partly suppressed in the 
FM/AO and AF/FO phases. In all these situations there is no possibility of joint 
spin-orbital fluctuations as the wave function of the GS is exactly known and has no 
quantum fluctuations in at least one of the two complementary subspaces. Under these 
circumstances separation of spin and orbital operators becomes exact and the GS is 
disentangled. This does not concern excited states [73], but in this section we are 
interested only in entanglement in the GS. Note, however, that at the special point 
x = y = — | one finds even three degenerate collective excitations: spin, orbital and 
spin-orbital wave, and the latter excitation is robust and does not decay into separate 
spin and orbital excitation [74 J. 

First we consider the variation of the full energy spectrum of a finite N = 4 site 
chain when the FM/FO state changes into the regime dominated by spin and orbital 
singlets. Therefere we study the SU(2)cg)SU(2) model ( 1251) along the symmetric line in 
the parameter space, p = x = y, see figure [2j Along this line interactions for S = \ 
spins and r = \ pseudospins appear on equal footing. At p = | one finds the high- 
symmetry SU(4) point and all three correlation functions are equal: ffl~9]) , ([20]) . 
and g ((Si ■ S i+ i)(Ti ■ T i+1 )) [75]. On the other hand, at p — | the model ( 125]) reads as a 
product of spin triplet and orbital triplet projection operators at each bond and its GS 
is exactly solvable — one finds two equivalent states with alternating spin and orbital 
singlets forming a spin-orbital valence bond (SOVB) phase [76J. These states may be 
obtained by a similar consideration as the Majumdar- Ghosh valence bond (VB) states 
in the ID frustrated J\ — J2 spin model with nearest J\ and next-nearest neighbour J2 
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Figure 3. Signatures of entangled spin and orbital operators in the SU(2)®SU(2) 
model (|25p. obtained by exact diagonalization of the ID chain of N = 8 sites at x = y 
with PBC: (a) the exact GS energy Eq (diamonds) and mean held energy Emf (squares) 
per bond; (b) intersite spin correlations Sij (|19p (circles), and composite spin-orbital 
correlations C,*j (j2Tj) (crosses); here orbital correlations (|20|) are degenerate with spin 
ones, Tij = S lJ . Entangled GS is found for p > -0.25 [76]. 

interactions, at J2 = \ J\ [ZZ|. The energy E = is given rigorously by alternating spin 
and orbital singlets along the chain. 

One finds a quantum phase transition (QPT) between the high spin-orbital FM/FO 
state (S t = It = 2) and the singlet entangled state (S t = % = 0) in the ID model ( l25l) . 
see figure [2l The QPT that occurs at p = — | is first order, as the energy levels shown 
in figure |2] cross and all intersite correlations change abruptly, see figure [3](b) . At this 
point the Hamiltonian (|25|) is a product of spin singlet and orbital singlet projection 
operators at each bond, so the FM/FO state has the lowest possible energy E = 0. 
But it suffices that either spin or orbital state is a triplet and thus the degeneracy of 
the GS is here much higher. The role of quantum fluctuations in the regime of p > | is 
easily recognized by considering further variation of the GS energy Eq with increasing 
p, see figure [2j Classically the energy of the model (125]) would be minimal at p = 0, but 
the quantum effects shift the energy minimum to the SU(4) point p — |; this state is 
nondegenerate. The energy Eq = is obtained again at p = |, where the Hamiltonian 
( |25|) is a product of spin triplet and orbital triplet projection operators at each bond - 
here one finds the SOVB phase explained above. 

While the GS energy E per bond of the ID spin-orbital model (l25l) is exactly 
reproduced by the MF energy normalized per one bond -Emf in the entire regime of 
singlet states (S t = % = 0) for p < — |, 

Emf = {(Si ■ £+1) + p) ((f, • f m ) + p) , (26) 

one finds large corrections to the MF energy for p > — |, see figure [3](a). This 
confirms that this MF decoupling procedure is not allowed as joint spin-orbital quantum 
fluctuations are then ignored which leads to uncontrollable errors in bond correlations. 

To characterize spin-orbital entanglement in the GS we evaluate intersite spin, 
orbital and joint spin-orbital bond correlations, defined in equations f lT§|) - fl2"l~|) . see 
figure Mjo)- Here we use s = | in equations ( !T9|) and (I2T|) . Note that spin and orbital 
correlations are intially weaker than the classical ones, i.e., = > — |, as they are 
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disturbed by the joint spin-orbital correlations ( 12 ip . The latter are negative in the 
entire regime of p > — \ and provide the dominating energy gain in the GS, including the 
SU(4) symmetric p = \ and the exactly solvable SOVB p = | points. It is of importance 
that spin-orbital entanglement is related to local properties of spins and orbitals on a 
bond. Therefore the entanglement phase diagram of a finite system is in agrement with 
the magnetic and orbital phase diagram of the infinite SU(2)cg>SU(2) model [75] . 

Another and a more precise quantity to quantify spin-orbital entanglement is von 
Neumann entropy defined as follows [79] . 

Sst = - Tr s {p 5 log 2 Ps} , (27) 

where ps = Tr T {\^)(^>\} is the reduced density matrix of the spin part in the state |\&) 
by integrating out all the orbital degrees of freedom by Tr r {- ■ •}. This measure captures 
well the correlation between the two types of degrees of freedom, and when spin and 
orbital operators factorize one finds Ss T = 0. It has been shown that von Neumann 
entropy makes a jump at the QPTs between the disentangled and entangled states [79] 
and may therefore be used to investigate the phase diagram of the SU(2)®SU(2) model. 
We suggest that QPTs in more complex systems of strongly correlated electrons may 
be investigated by calculating von Neumann entropy in the future. 

3.3. Entanglement in tig spin-orbital models 

Entanglement plays an important role in realistic spin-orbital models for t 2 g orbital 
degrees of freedom which may be considered as being more quantum than the respective 
models for e g electrons. This is a consequence of two orbital flavours being active along 
each cubic axis, while in e g systems there is de facto just one directional orbital which is 
active along a given cubic axis while the other orbital is inactive, so these latter models 
are more classical. 

The best known example of transition metal oxides with the physical properties 
controlled by spin-orbital entanglement are the -RVO3 perovskites. Two magnetic phases 
compete with each other at low temperature, and one finds the C-AF phase accompanied 
by C-AO order in compounds with a large ionic radius tr of R ions, i.e., for Pr,- • -,La, 
while for R ions with smaller ionic radius the G-AF phase accompanied by C-AO order is 
more stable at T = and the C-AF phase occurs only in a window of finite temperature. 
In these phases both the magnetic moments and the occupied orbitals alternate in the 
oh planes along both cubic axes, but the order along the c axis is different, as shown in 
figure HI Note that along the c axis the Goodeneough-Kanamori rules are obeyed. The 
situation looks different for the ah planes, where the magnetic moments and occupied 
orbitals alternate, see also section 6, but one should keep in mind that inactive c orbitals 
are singly occupied at every site and rather strong AF superexchange arises along both 
the a and b axis due to their excitations. Thus, one may classify this FO order of 

c orbitals accompanied by AF order of S = 1 spins, and once again the Goodeneough- 
Kanamori rules are followed. 
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Figure 4. Schematic view of the Goodenough-Kanamori rules on the example of two 
AF phases observed in the i?VO*3 perovskites, with complementary magnetic (arrows) 
and orbital (squares) order of active t2 g orbitals {a, b} in the ac plane: (a) C-AF spin 
order accompanied by G-AF orbital order; (b) G-AF spin order accompanied by G-AF 
orbital order. Both spins and orbitals alternate along the b axis (not shown). This 
figure is reproduced from 



Although the above discussion shows that the GS of the -RVO3 compounds does not 
include spin-orbital entanglement, it is of interest to investigate the spin orbital models 
for the perovskite systems with active t2 g orbitals: (i) the titanate model (d 1 model) 
valid for the i?Ti03 perovskites [60JET], and (ii) the vanadate model (d 2 model) which 
describes the _RV0 3 perovskites [67J[7T]. Spin-orbital entanglement arises along the c 
axis, where both active t2 9 orbitals {a, b} contribute and may lead to entangled states. 
To avoid additional complications due to partly occupied orbitals, one may assume that 
the c orbitals are empty at every site in the d 1 model, while they are occupied in the d 2 
model — in both cases they cannot lead to any entangled states. 

The spin-orbital models for the cubic perovskites with active t<i g orbital degrees 
of freedom at either Ti 3+ or V 3+ ions are of the general form given in equation ( 1131) . 
The orbital operators and are rather complex and depend on the multiplet 
structure of the Ti 2+ and V 2+ excited states, respectively. They include the terms which 
break the SU(2) symmetry of the orbital interactions present at 77 = 0, reducing it to the 
cubic symmetry. Their explicit form may be found in the original publications. Here 
we give only the simplified SU(2)<S>SU(2) form of the interactions for the bonds along 
the c axis, 

Hj = \ J H ■ + ^) ■ n+i + I) , (28) 

i 

where spin-orbital entanglement is expected. For S — \ and S = 1 the above general 
form reproduces the respective limits obtained for the d 1 and d 2 models at 77 = 0; 



Fingerprints of spin-orbital entanglement in transition metal oxides 



19 



otherwise it is approximate. We assume below that 



n ia + n ib = 1 



(29) 



for t-ig orbitals (jHJ) in case of the d 1 model, and ni C = 0. 

In the .RVO3 perovskites, the crystal-field splitting breaks the cubic symmetry in 
distorted V0 6 octahedra, as suggested by the electronic structure calculations |80j and 
derived using the point charge model [81]. One finds again the same filling of {a, b} 
orbitals as given in equation ( l2"9~j) . but n ic = 1. This defines the t 2g orbital degrees 
of freedom in both cases as {a, b} = {yz, zx} orbitals along every cubic direction, 
and the superexchange (|40l) contains the orbital operators fj = {r^, rf, t[} (and their 
components). 

A method of choice to demonstrate spin-orbital entanglement is here again exact 
diagonalization of finite chains, performed for both the d 1 and d 2 model pQ. In the 
d 1 model the Hamiltonian at 77 = reduces to the SU(4) model, and indeed all three 
bond correlation functions are equal for N = 4 sites [1J, S^ = = = —0.25 as 
shown in figure [5(a). For larger systems these correlations are also equal but somewhat 
weaker and one finds S^ = TJy = Cy = —0.21502 in the thermodynamic limit [75J. 
By a closer inspection one obtains that the GS wave function for the four-site cluster is 
close to a total spin-orbital singlet, involving a linear combination of (spin singlet/orbital 
triplet) and (spin triplet/orbital singlet) states for each bond (ij). This result manifestly 
contradicts the celebrated Goodenough-Kanamori rules [14J, as both spin and orbital 
correlations have the same sign. When rj increases, the charge fluctuations which 
contribute to the superexchange concern different states in the multiplet structure breaks 
the SU(4) symmetry — one finds that the bond correlations are then different and 
< Cij < < in the regime of spin singlet (S t = 0) GS. Here the Goodenough- 
Kanamori rule which suggests complementary spin/orbital correlations is still violated. 

The vanadate d 2 model (for 5=1 spins) [71] behaves also in a similar way in a 
range of small values of 77, with all three S^, and correlations being negative, 
see figure [5(b). Most importantly, the composite spin-orbital correlations are here finite 
(Cij < 0) which implies that spin and orbital variables are entangled, and the MF 
factorization of the GS into spin and orbital part fails. In this regime the spin and orbital 
correlations are both negative and contradict the Goodenough-Kanamori rules [14J of 
their complementary behaviour. Only for sufficiently large 77 do the spins reorient at 
the QPT to the FM GS, and decouple from the orbitals. In this regime, corresponding 
to the experimentally observed C-AF phase of LaVOs (and other cubic vanadates), 
spin-orbital entanglement ceases to exist in the GS. However, as we will see below, it 
has still remarkable consequences in experiments at finite temperature, where entangled 
spin-orbital excited states contribute and decide about the thermodynamic properties. 

A crucial observation concerning the applicability of the Goodenough-Kanamori 
rules to the quantum models of t 2g electrons in one dimension can be made by comparing 
spin exchange constants Jy calculated from the exact GS |<3>), 



J« = <*l4 7) |*>, 



(30) 
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Figure 5. Evolution of: (a-b) intersite spin Sij (fl9|) (filled circles), orbital Ty 
(|20| (empty circles), and spin-orbital Cy (|2"Tj) (x) correlations, and (c-d) exchange 
constants J t j (|14[) along the c axis, obtained by exact diagonalizaton of spin-orbital 
model on a ID chain of N — 4 sites with periodic boundary conditions, for increasing 
Hund's exchange i]. Panels (a,c) for S = 1/2; panels (b,d) for 5 = 1. In the shaded 
areas of (c) and (d) the negative spin correlations < do not follow the sign of the 
exchange constant < 0, and the classical Goodenough-Kanamori rules are violated. 
This figure is reproduced from [llj . 



with intersite spin correlations ( |T9|) obtained also exactly. One finds that exchange 
interaction which is formally FM ( < 0) in the orbital-disordered phase at low values 
of r) [see figure [5]Jc) and figure EJ^d)] is in fact accompanied by AF spin correlations 
(Sij < 0), so JijSij > and the GS energy calculated in the MF theory is de facto 
enhanced by this term PQ, contrary to what happens in reality. 

In contrast, similar analysis (not shown) performed for the d 9 spin-orbital model 
(jl~3]) derived for Cu 2+ ions with e g orbital degrees of freedom in KCuF 3 [70], gave 
correctly JijSij < 0. Hence, in spite of enhanced quantum fluctuations of the spin- 
orbital nature [22] , one finds here that spin correlations follow the sign of the exchange 
constant [I]. This remarkable difference between ti g and e g systems originates from 
composite spin-orbital fluctuations, which are responsible for the 'dynamical' nature 
of exchange constants in the former case. They exhibit large fluctuations around the 
average value, measured by the variance, 

u = {m(j^Y\$) - ji) 112 . (si) 

Here again the average is calculated exactly from the Lanczos diagonalization of a finite 
chain of length N = 4 sites. 

As an illustrative example, we give here the values found in the d l and d 2 model at 
rj — pQ. While the average spin exchange constant is small in both cases ( ~ for d 1 , 
Jij ~ —0.04 for d 2 ), fluctuates widely over both positive and negative values. In the 
d 1 model the fluctuations between (S = 0/T = 1) and (S = 1/T = 0) wave functions on 
the bond are so large that 5 J = 1 ! They survive even quite far from the high-symmetry 
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SU(4) point (at 77 ~ 0.1), and stabilize spin-orbital singlet phase in a broad regime of 
7]. Also in the d 2 model the orbital bond correlations change dynamically from singlet 
to triplet, resulting in 5 J > \J lj \, with 5 J = \{l - (2T y + \) 2 } 1/2 ~ 0.25, while these 
fluctuations are small for d 9 model involving e g orbitals, see also section 15.11 

We emphasize that spin and orbital correlations on the bonds, as well as composite 
spin-orbital correlations which occur in spin-orbital entangled states for realistic 
parameters with finite Hund's exchange, determine the magnetic and optical properties 
of titanates and vanadates. These correlations change with increasing temperature as 
then also excited states contribute and decide about their thermal evolution. Therefore, 
the correct theoretical description of experimental results is challenging and requires 
adequate treatment of excited states which captures their essential features, including 
their possible entanglement. This makes simple approaches based on the MF decoupling 
of spin and orbital operators not reliable and requires either exact diagonalization 
of finite clusters or advanced numerical methods such as multiscale entanglement 
renormalization ansatz (MERA) [82]. 

In the next section we show that composite spin-orbital fluctuations play a crucial 
role in the .RVO3 perovskites. They are responsible for the temperature dependence of 
the optical spectral weights in LaV0 3 [35], contribute to the remarkable phase diagram 
of these systems [81] and trigger spin-orbital dimerization in the C-AF phase of YVO3 
in the intermediate temperature regime [83] . Remarkably, all these properties including 
the observed dimerization in the magnetic excitations may be seen as signatures of spin- 
orbital entanglement in the excited states which becomes relevant at finite temperature. 

4. Entangled states in the RVO3 perovskites 

4-1. Optical spectral weights for LaVO^ 

The coupling between spin and orbital operators in the spin-orbital superexchange may 
be so strong in some cases that it leads to a phase transition modifies the magnetic 
order and excitations at finite temperature — an excellent example of this behaviour are 
the RVO3 perovskites, as explained below. Although the C-AF phase, observed in the 
entire family of KVO3 compounds [21], where R=Lu,- ■ -,La stands for a rare earth atom, 
satisfies to some extent the Goodenough-Kanamori rules [H], with FM order along the c 
axis where the active a and b orbitals dHJ alternate — the G-AO order of {a, b} orbitals 
is very weak here and the orbital fluctuations play a very important role [67]. This 
situation is opposite to the frozen and classical AO order in LaMnOs, which can explain 
both the observed magnetic exchange constants [13] and the distribution of the optical 
spectral weights [H]. In LaV03 the FM exchange interaction is enhanced far beyond 
the usual mechanism following from the splitting between the high-spin and low-spin 
states due to finite Hund's exchange Jh- Evidence of orbital fluctuations in the -RVO3 
perovskites was also found in pressure experiments, which show a distinct competition 
between the C-AF and G-AF spin order, accompanied by the complementary G-AO 
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and C-AO order of {a, b} orbitals [2Tj . 

The spin and orbital order along the c axis are not entangled in the GS of the 
-RVO3 perovskites (due to either FM or FO order), but entangled states contribute 
at finite temperature. As the first manifestation of spin-orbital entanglement in the 
cubic vanadates at finite temperature we discuss briefly the evaluation of the low- 
energy optical spectral weight from the spin-orbital superexchange for LaVOs, following 
equation (ITS]). The superexchange operator 7-Lj (fTT|) is here considered for a bond 
(ij) || 7, and arises as a superposition of individual dfd 2 ^ d\d^ charge excitations to 
different spin states in the upper Hubbard subbands labeled by n |35j. One finds the 
superexchange terms H^- for a bond (ij) along the c axis, 

Hi% = - \jt x (2+3.3) Q - n ■ , (32) 

Hi% = ~ (l - S,rS 3 ) Q v,r: - r;+;-.o7:'7f) , (33) 

= - \jn (l -Si-ty Q + rfr; - rfrj) , (34) 

and if^- for a bond in the ah plane, 

H n£ = ~ \jri(2 + S t -S 3 ) Q - rlr^ , (35) 

B &=-h J ( 1 ~ 3t 'fy {j TT * TT >- \ T ^) ' (36) 

H Sh - tq Jt * i 1 - Q T r ^ r ^ + 2r ^) • (37) 

When the spectral weight is evaluated following equation ffTB"]) . it is reasonable to try first 
the MF approximation and to separate spin and orbital correlations from each other. 
The spectral weights require then the knowledge of spin correlations (fT9l) : along the 
c axis in (j3"2^ - (j3"4"j) . and within the ah planes in (|3"5|) - (|3"7|) . as well as the corresponding 
intersite orbital correlations (fj ■ fj) and (rfr"), with a = x,y, z. From the form of the 
above superexchange contributions one sees that high-spin excitations H^J- support the 
FM coupling while the low spin ones, and H^, contribute with AF couplings. 

The high-spin spectral weight (Tl8|) in the MF approximation is given by 

< F = pn^-Sj +2} (± -n ■ . (38) 

The low-energy optical spectral weight for the polarization along the c axis K[°l xp 
decreases by a factor close to two when the temperature increases from T ~ to T = 300 
K [25] - this change is much larger than the one observed in LaMnOa [16]. However, 
the theory based on the MF decoupling of the spin and orbital degrees of freedom gives 
only a much smaller reduction of the weight close to 33% when a frozen orbital order 
(similar to LaMnOa) with (fi-fj) = — | is assumed and has to fail in explaining the 
experimental data [T5] . In spite of weakening spin and orbital intersite correlations when 
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Figure 6. Optical spectral weights (fT5|) for the high-spin optical subband (n = 1; 
solid line), low-spin optical subbands (n = 2,3; dashed lines) and total (dotted 
line) obtained from the excitations to different Hubbard subbands along the c axis 
([3"2"]) - (fM)l . Filled circles show the effective carrier number (in the energy range 

oj < 3 eV which corresponds to the high-spin excitations) for LaV03, presented in Fig. 
5 of [25]. Dashed-dotted line shows the spectral weight K± obtained from the MF 
decoupling <j3Sj). Parameters: rj = 0.f2, V c = 0.9J, V ab = 0.2J. 



the temperature increases in the experimental range < T < 300 K, this variation is 
clearly not sufficient to describe the experimental data. Instead, when both spin and 
orbital correlations are reduced and vanish above the common transition temperature 
T N1 = T o, the MF theory predicts that the spectral weight decrease fast and do not 
change above Tni, see figure EJ contrary to experiment. 

On the other hand, when a cluster method is used to determine the optical spectral 
weight from the high-spin superexchange term (|32|) by including orbital as well as 
joint spin-and-orbital fluctuations along the c axis, the temperature dependence of the 
spectral weight resulting from the theory persists above Tni and follows the experimental 
data [35], see figure El In this approach a cluster of N = 4 sites is solved exactly with the 
MFs originating from the spin and orbital order below Tni and Too, respectively, and 
a free cluster is solved at high temperature when the long-range order vanishes. This 
reflects the realistic temperature dependence with spin and orbital correlations being 
finite in this latter regime, in contrast to the one-site MF approach. 

The satisfactory description of the experimental data shown in figure [6] may be 
considered as a remarkable success of the theory based on the spin-orbital superexchange 
model derived for the RVO3 perovskites. It proves that spin-orbital entangled states 
contribute in a crucial way in the entire regime of finite temperature. In addition, the 
theoretical calculation predicts that the high energy spectral weight (n — 2, 3) is low for 
the polarization along the c axis. The spectral weight in the ab planes behaves in the 
opposite way — it is small at low energy, and large at high energy (but not as large as 
the low-energy one for the c axis). This weight distribution and its anisotropy between 
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the c and ab directions reflects the nature of spin correlations on the bonds, which are 
FM and AF in these two directions. We expect that future experiments will confirm 
these theory predictions for the ab polarization. 

4-2. Phase diagram of the RVO3 perovskites 

The phase diagram of the i?Mn03 perovskites [84] indicates that spin and orbital 
energy scales separate which makes it possible to describe the experimental data for the 
magnetic exchange constants and the optical spectral weights using the disentangled 
spin-orbital superexchange [13, 16J. In contrast, in the .RVO3 perovskites the phase 
diagram suggests the proximity of spin and orbital energy scales [21J. Experimental 
studies have shown that the C-AF order is common to the entire family of the 
-RVO3 vanadates, and in general the magnetic transition occurs soon below the orbital 
transition when the temperature is lowered further, i.e., Tni < Too- LaVC>3 is an 
exception here and these transitions occur almost simultaneously, with Tni — Too |21j. 

When the ionic radius r R decreases, the Neel temperature T N1 also decreases, while 
the orbital transition temperature Too increases, passes through a maximum close to 
YVO3, and next decreases when LuV0 3 is approached [21]. One finds that the C-AF 
order develops in LaV03 below Tvi — 143 K, and is almost immediately followed by a 
weak structural transition stabilizing the weak C-AO order at Too — 141 K [21]. This 
provides a constraint on the theoretical model and is an experimental challenge to the 
theory which was addressed using the spin-orbital superexchange model extended by 
the orbital- lattice coupling [81]. 

In order to unravel the physical mechanism responsible for the surprising decrease 
of Tqo from YVO3 to LUVO3 one has to analyse in more detail the evolution of GdFeC>3 
distortions for decreasing ionic radius r# [ST]. Such distortions are common for the 
perovskites [85], and one expects that they should increase when the ionic radius 
r R decreases, as observed in the RMn0 3 perovskites [SIT . In the -RVO3 family the 
distortions are described by two subsequent rotations of VOg octahedra: (z) by an angle 
i? around the b axis, and (ii) by an angle ip around the c axis. Increasing angle d causes 
a decrease of V-O-V bond angle along the c direction, being ir — 2a9, and leads to an 
orthorhombic lattice distortion u = (b — a) /a, where a and b are the lattice parameters 
of the Pbnm structure of -RVO3. By the analysis of the structural data for the .RVO3 
perovskites [S61IS7] one finds the following empirical relation between the ionic radius 
tr and the angle d: 

tr = r — a sin 2 $ , (39) 

where ro = 1.5 A and a = 0.95 A are the empirical parameters. This allows one to 
use the angle to parametrize the dependence of the microscopic parameters of the 
Hamiltonian and to investigate the transition temperatures Too and Tni as functions 
of the ionic radius r R . 

The spin-orbital model introduced in [81] to describe the phase diagram of RVO3 



Fingerprints of spin-orbital entanglement in transition metal oxides 



25 



was thus extended with respect to its original form [57] and reads: 

U = J E { (4 7) ^4 + S 2 ) + + E^^e^Tl 

(y>ll7 * 



Y^t* + \nK{u-u q {$)}\ (40) 



where 7 = a, 6, c labels the cubic axes, and the orbital operators and /Qj are given 
in [71]. The superexchange oc J is supplemented by the crystal field term oc E z , the 
orbital interaction terms oc V c and oc V a b induced by lattice distortions, and the orbital- 
lattice term oc g which is counteracted by the lattice elastic energy oc K. All these terms 
are necessary in a realistic model [HI] to reproduce the complex dependence of the orbital 
and magnetic transition temperature on the ionic size in the -RVO3 perovskites. 

The crystal-field splitting oc E z between a and b orbital energies in equation (liOjl 
is given by the pseudospin rf operators, 

r- = ^(n ia - n ib ) , (41) 

which refer to two active orbital flavors {a, b} in RVO3. It is characterized by the 
vector Q = (tt, tt, 0) in reciprocal space and favours the C-AO order. Thus, this 
splitting competes with the (weak) G-AO order supporting the observed C-AF phase 
at temperature T < T$i, effectively weakening this type of magnetic order. 

As for instance in LaMnOs, the orbital order in the -RVO3 perovskites arises 
due to joint orbital interactions which are a superposition of the superexchange and 
interactions induced by lattice distortions. These latter terms are twofold: (i) intersite 
orbital interactions {V a b, V c } (which originate from the coupling to the lattice), and (ii) 
orbital-lattice coupling oc g which induces orbital polarization (rf ) 7^ for finite lattice 
distortion u. The orbital interactions induced by the distortions of the VC>6 octahedra 
and by the GdFe0 3 distortions of the lattice, V a t, > and V c > 0, also favour the C-AO 
order (like the crystal field term with E z > 0). Note that V c > counteracts the orbital 
interactions included in the superexchange via operators. 

The last two terms in equation fj4*0|) are particularly important for ions with small 
ionic radii Tr. They describe the linear coupling oc g > between active {a, b} orbitals 
and the orthorhombic lattice distortion u. The elastic energy which counteracts lattice 
distortion u is given by the force constant K, and iV is the number of V 3+ ions. The 
coupling 

^cfr = gu (42) 

may be seen as a transverse field in the pseudospin space which competes with the Jahn- 
Teller terms {V a b, V c }. While the eigenstates of rf operator, ± \b)), cannot be 

realized due to the competition with all the other terms, increasing lattice distortion u 
(increasing angle d) gradually modifies the orbital order and intersite orbital correlations 
towards this type order. 



Fingerprints of spin-orbital entanglement in transition metal oxides 



26 



Except for the superexchange parameter J, all the parameters in the extended 
spin-orbital model fj4*0]) depend on the tilting angle In case of V c one may argue 
that its dependence on the angle ■d is weak, and a constant = 0.26J was chosen 

in [81] in order to satisfy the experimental constraint that the magnetic and orbital order 
appear almost simultaneously in LaV03 [21] . The experimental value T^ p = 143 K for 
LaV0 3 [21] was fairly well reproduced in the present model taking J = 200 K, the value 
which is also consistent with the magnon energy scale [26]. The functional dependence 
of the remaining two parameters {E z ($), Vabi'd)} on the tilting angle $ was derived from 
the point charge model [81] using the structural data for the i?V0 3 series [ES1IEI] — one 
finds: 

E z (ti) = Jv z sin 3 ?? cost?, (43) 

Vabiti) = J v ab sin 3 cos $ . (44) 

Finally, the effective coupling to the lattice distortion g e E{$) (H2]) has to increase faster 
with the ncreasing angle $ as otherwise the nonmonotonous dependence of Too on 
$ (or on the ionic radius r#) cannot be reproduced by the model, and the following 
dependence was shown [Hi] to give a satisfactory description of the phase diagram of 
the -RVO3 perovskites: 

9efi($) = Jv g sin 5 $ cost? . (45) 

Altogether, magnetic and orbital correlations described by the spin-orbital model 
( f40l) . and the magnetic Tni and orbital Too transition temperatures, depend on three 
parameters: {v z , v a b, v g }. 

Due to the spin-orbital entanglement which is activated by finite temperature in the 
-RVO3 perovskites, it is crucial to design the MF approach in such a way that the spin- 
orbital coupling is described beyond the factorization of spin and orbital operators. As 
usually, the correct MF treatment of the orbital and magnetic phase transitions in the 
-RVO3 perovskites requires the coupling between the on-site orbital, (t 2 ) g = ||(t* — t?)|, 
and spin order parameters in the C-AF phase, {S-)c, as well as including a composite 
spin-orbital (Sfr?) order parameter, similar to that introduced before for the i?Mn0 3 
perovskites [12]. However, the on-site MF approach including the above coupling [HE] 
does not suffice for the -RVO3 compounds as the orbital singlet correlations (Ti-Tj) on the 
bonds along the c axis play here a crucial role in stabilizing the C-AF phase [67J and the 
orbital fluctuations are also important [89J . Therefore, the minimal physically acceptable 
approach to the present problem is a self-consistent calculation od spin and orbital 
correlations for an embedded bond (ij) along the c axis, coupled by the MF terms to its 
neighbours along all three cubic axes [HI] • This procedure, with properly selected model 
parameters, led to the successful description of the experimental phase diagram [2TJ , see 
figure [71 The evolution of the orbital correlations with varying temperature and with 
decreasing r# plays a prominent role in the success of the theoretical description of the 
phase diagram of the -RVO3 perovskites. One finds that indeed the orbital order occurs 
first at a higher temperature and is followed by the magnetic transition in the i?V0 3 
perovskites with a smaller ionic radius r^, to the left of LaVOs. 



Fingerprints of spin-orbital entanglement in transition metal oxides 27 

Lu Y Sm Pr La 

250 r— — . ■ ■ — — — — — i 




' — 1 — — — — — — 1 

1.1 1.2 1.3 1.4 

r R (A) 



Figure 7. The orbital transition temperature Too anc ^ Neel transition temperature 
Tni for the onset C-AF order as functions of varying ionic size in the -RVO3 
perovskites, as obtained in experiment (full and empty circles) and from the theory 
(solid lines). Dashed lines show Too and Tni obtained under neglect of orbital-lattice 
coupling (at g = 0). The inset shows the GdFeC>3-type distortion, with the rotation 
angles ■& and ip corresponding to the data of YVO3 [86j . This figure is reproduced 
from [81]. 

The non-monotonous dependence of the orbital transition temperature Too 011 
the ionic radius r# may be understood as follows. Too increases first with decreasing 
ionic radius as the Jahn- Teller term in the ab planes, V^,(??), increases and induces 
the orbital correlations which stabilize the G-AO order. The coupling g e s{^) to the 
lattice ( 145|) is then rather weak, with (rf ) ~ 0.03 in LaV0 3 , but g e $ increases faster 
than the interaction V^($) ( 1441) . Finally, the former term dominates and the G- 
AO order parameter is almost equal to the competing with it "transverse" moments, 
{rf) — (t[)g- Therefore, the G-AO order gets weaker and the transition temperature 
Too is reduced. Note that in the entire parameter range the orbital order parameter 
{rf) ~ {^1)0 is substantially reduced from the classical value {r^cm&x — \ by singlet 
orbital fluctuations, being for instance {rf)a — 0.32 and 0.36 for LaV0 3 and LuV0 3 . 

It is quite remarkable that the magnetic exchange constants {J a b, J c } are modified 
solely by the changes in the orbital correlations described above. The superexchange 
constant J does not change and the reductions of Tni with decreasing follows only 
from the evolution of the orbital state [81J. One finds that also the width of the magnon 
band, given byWc-AF — 4(J a & + | J c |) at T = 0, is reduced by a factor close to 1.8 from 
LaV0 3 to YV0 3 . This also agrees qualitatively with surprisingly low magnon energies 
observed in the C-AF phase of YV0 3 [26J. 

Summarizing, the microscopic model (HO]) is remarkably successful in describing 
gradual changes of the orbital and magnetic correlations under increasing Jahn- Teller 
interactions and the coupling to the lattice which both suppress the orbital fluctuations 
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along the c axis, responsible for rather strong FM spin-orbital superexchange [67J. 
It describes well the systematic experimental trends for both orbital and magnetic 
transitions in the -RVO3 perovskites |8T] , and is able to reproduce the observed non- 
monotonic variation of the orbital transition temperature Too f° r decreasing ionic radius 
tr. Another consequence of the spin-orbital entanglement in the perovskite vanadates, 
the spin-orbital dimerization along the c axis in YVO3, is shortly discussed in the next 
subsection. 

4-3. Peierls dimerization in YVO3 

The third and final example of the spin-orbital entanglement at finite temperature in 
the family of vanadate perovskites is the existence of a remarkable first order magnetic 
transition at Tn 2 = 77 K from the G-AF to the C-AF spin order with rather exotic 
magnetic properties, found in YVO3 [90J. This magnetic transition is surprising and 
rather unusual as the staggered moments are approximately parallel to the c axis in the 
G-AF phase, and reorient above Tn 2 to the ab planes in the C-AF phase, with some small 
alternating G-AF component along the c axis. First, while the orientations of spins in 
C-AF and G-AF phase are consistent with the expected anisotropy due to spin-orbit 
coupling [83], the observed magnetization reversal with the weak FM component remains 
puzzling and found no explanation in the theory so far. Second, it was also established 
by neutron scattering experiments [26] that the energy scale of magnetic excitations is 
considerably reduced for the C-AF phase (by a factor close to two) as compared with 
the magnon dispersion measured in the G-AF phase. The magnetic order parameter 
in the C-AF phase of LaVOs is also strongly reduced to ~ 1.3/xb, which cannot be 
explained by rather small quantum fluctuations in the C-AF phase [91]. Finally, the 
C-AF phase of YVO3 is dimerized. Until now, only this last feature found a satisfactory 
explanation in the theory [92], see below. 

The observed dimerization in the magnon spectra in YVO3 motivated the search for 
its mechanism within the spin-orbital superexchange model. Dimerization of AF spin 
chains coupled to phonons is well known and occurs in several systems [93] . In the spin- 
orbital model for the -RVO3 perovskites a similar instability might also occur without 
the coupling to the lattice when Hund's exchange is sufficiently small. In particular, 
the GS at rj = may be approximated by the dimerized chain with strong FM bonds 
alternating with the AF ones, if such chains are coupled by AF interactions along the a 
and b axes [94J (the ID chain would give then the entangled disordered GS). For finite 
and realistic rj ~ 0.13 the chain is FM (due to the weak coupling to the neighbouring 
chains in ab planes) [71] and at first instance any dimerization appears surprising. 

Before addressing the question of magnon excitations in the C-AF phase of 
YVO3 stable at intermediate temperature, let us consider first the ID spin-orbital 
superexchange model along the c axis, as in YVO3. The Hamiltonian is given by [92] , 






(46) 
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where Sj represent S = 1 spins and Tj are r — h orbital pseudospins, respectively, and 
is a constant proportional to Hund's exchange which stabilizes FM spin correlations. 
This expression is somewhat simplified with respect to the full spin-orbital model for 
YV0 3 [71], but reflects its essential features and guarantees that the GS is FM when 
7# ~ 0.1. The FM GS state is disentangled — it is allowed to use the MF decoupling pQ, 
and to decompose the above Hamiltonian (J36]) into the spin (Hs) and orbital (H T ) part, 
Hsr — Hs + H T . This disentangled chain may be now studied either by density-matrix 
renormalization group applied to transfer matrices (TMRG) [95], or by an analytical 
approach, the so-called modified spin-wave theory of Takahashi [96J. 

It is easy to understand why the spin-orbital dimerization occurs at finite 
temperature. The crucial concept is the interrelation between spin and orbital 
correlations in the ID spin-orbital chain: spin correlations determine the exchange 
interactions in the orbital channel H T , and the orbital ones are responsible for the 
spin exchange in Hs- In the GS the spin state is rigid, and spin correlations on the 
bonds + 1) along the c axis are saturated, i.e., (Sj ■ Sj+i) = 1, and do not allow 
for any alternation in the orbital interactions which are determined by them. But when 
temperature increases the thermal fluctuations soften the FM order and the spin-orbital 
chain may dimerise [92]. Important here is the rather dense spectrum of low energy 
excited states in the spin-orbital chain [97], which are entangled and all contribute to 
the thermal averages used to calculate spin and orbital correlations. We emphasize that 
the dimerization in the spin-orbital chain may be seen as a signature of entanglement in 
excited states in the C-AF phase which contribute at finite temperature. The exchange 
constants alternate along the c direction between a stronger {Jc 1 ^) and weaker (J7"i 1 ' ) ) 
exchange with 5s > 0, 

jQ)=j e (l + 6 s ), jj?)=J c (l-5s). (47) 

Similar expressions are also found for the orbital exchange interactions which favour AO 
order and have alternating strength with S T > 0, 

jW=J T (l + S T ), jW = j T (i-5 T ). (48) 

While the spin and orbital operators are disentangled in the FM ground state, 
one may consider a coupled FM spin chain to an orbital chain with interactions which 
favour weak AO order accompanied by orbital fluctuations, as realized in the C-AF 
phase. The spin (orbital) exchange interaction along the chain is then determined by 
the bond orbital (spin) correlations. They are defined as follows: 

Jc = ^(r r T j+1 + Tj -T^-l) + ~ - 7ff , (49) 

J T = l(Sj-S j+1 + Sj-Sj-x) + 1 , (50) 

and have to be determined self-consistently, together with spin and orbital correlations 
along the chain. For the parameters selected in equation ( |46l) . one finds that FM spin 
correlations with J c < are accompanied by J T > that favours AO order along the 
chain. Such complementary spin and orbital correlations are indeed found in the entire 
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Figure 8. (Color online) Phase diagram with two undimerized (U) phases and a 
dimerized (D) phase in between, and dimerization parameters 6s and S T , obtained for 
the spin-orbital model in the MF decoupling. Dotted vertical lines indicate the 
onset of a dimcriscd phase under increasing/decreasing temperature. For more details 
see [22]) where these data were presented and discussed. 



temperature regime, in agreement with the Goodenough-Kanamori rules. But in the 
intermediate temperature range, when the spins start to fluctuate and their correlations 
are not rigid anymore, the dimerization sets in, see figure El Finite temperature is here 
essential as dimerized spin correlations support then the dimerized orbital correlations. 
Hence, the dimerization occurs here simultaneously in both channels and has a dome- 
shaped form, with a maximum at T ~ 0.2 J [92]. The dimerization in the FM chain is 
much stronger than the one in the AO chain, but they have to coexist in the present 
self-consistent treatment of this phenomenon. When temperature is high enough, in the 
present case for T > 0.49 J, the dimerization vanishes again as the spins and the orbitals 
are disordered by thermal fluctuations. The phase transition at finite temperature 
between a uniform and a dimerised phase, shown in figure [BJ is a consequence of the 
MF decoupling. 

The microscopic model PS"]) , which explains the anisotropy in the exchange 
constants ( )47l) as following from the joint dimerization that occurs in the spin-orbital 
chain with FM spin order at finite temperature |92j, helps to understand the magnon 
dispersion found in YVO3 by the neutron scattering [26]. The observed spin- wave 
dispersion may be explained by the following effective spin Hamiltonian for the C-AF 
phase, derived assuming again that the spin and orbital operators may be disentangled 
which is strictly valid only at T = 0: 

H S = J C {! + (-!)%} Si -S i+1 

{i,i+l)\\c 

+ J ab St-Sj + K^iS*) 2 . (51) 

(ij)\\ab i 
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Figure 9. (Color online) Spin-wave dispersions (solid lines) as obtained in the 
LSW theory along the representative directions in the Brillouin zone for the dimerized 
C-AF phase of YVO3 with experimental exchange constants [26 : J a b = 2.6 meV, 
J c = 3.1(1 ± 6s) meV and 5s = 0.35. Parameters: J = 30 meV, 77 = 0.13, 5s — 0.35, 
K z = 0. The experimental points |26j measured by neutron scattering at T — 85 K 
are reproduced by circles. The high symmetry points are: T = (0,0,0), M = (jt, it, 0), 
R = (it, 7r, 71"), Z = (0, 0, 7r). This figure is reproduced from [71] . 




Following the linear spin-wave (LSW) theory the magnon dispersion is given by 



with 



u ± (k) = 2J^2J ab +\J c \ + ^K z ±J cV l/ 2 ^ -(2J afe7k ) 2 , (52) 



7k = - (cos k x + cos k y ) , (53) 



2 

77k = cos 2 k z + 6g sin 2 k z . (54) 

The single-ion anisotropy term oc K z is responsible for the gap which opens in spin 
excitations. Two modes measured by the neutron scattering [26J are well reproduced by 
0J±(k) given by equation ( 152]) when the experimental exchange constants are inserted: 
J a b = 2.6 meV, J c = —3.1 meV, 5s = 0.35, see figured (further improvement including 
a finite gap at the F point are obtained taking finite K z > [26J. This shows that 
while the essential features seen in the experiment are well reproduced already by the 
present simplified spin exchange model (EE]), the spin interactions are more complex 
in reality [25] . 

Summarizing, spin-orbital entanglement in the excited states is also responsible 
for the exotic magnetic properties of the C-AF phase of YVO3. They arise from the 
coupling between the spin and orbital operators which triggers the dimerization of the 
FM interactions ( |47j) as a manifestation of a universal instability which occurs in FM 
spin chains at finite temperature, either by the coupling to the lattice or to purely 
electronic degrees of freedom [92J. This latter mechanism could play a role in many 
transition metal oxides with (nearly) degenerate orbital states. 
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5. Entanglement in the ground states of spin-orbital models 

5.1. Kugel-Khomskii model 

As shown in section 3, the GSs of certain spin-orbital models are entangled and this 
will likely influence future experimental studies. As an example we discuss here the 
Kugel-Khomskii d 9 model on a bilayer and analyse the d 1 spin-orbital on a frustrated 
triangular lattice in section 15.21 While the coexistence the A-type AF (A-AF) order and 
the C-AO order below Tn — 39 K is well established in the KCuF 3 perovskite [981199] 
and this phase is well reproduced by the spin-orbital superexchange d 9 model [70], the 
model poses an interesting theoretical question: which types of coexisting spin-orbital 
order (or disorder) are possible when its microscopic parameters are varied? So far, 
it was only established that the long-range AF order is destroyed by strong quantum 
fluctuations [23l|23J[70] , and it has been shown that certain spin disordered phases with 
VB correlations may be stabilized by local orbital correlations [22,100j. However, the 
phase diagram of the Kugel-Khomskii d 9 model was not studied systematically beyond 
the MF approximation or certain simple variational wave functions and it remains an 
outstanding problem in the theory [22] . 

The simplest spin-orbital models are obtained when transition metal ions are 
occupied by either one electron (m = 1), or by nine electrons (m = 9); in these 
cases the Coulomb interactions ([9]) contribute only in the excited states (in the 
d 2 or the d 8 configuration) after a charge excitation between two neighboring ions, 
d™d™ ^ d T [ l+ d™~ . A paradigmatic example of the spin-orbital physics is obtained in 
the case of a single hole in the d shell, as realized for the d 9 (m = 9) configuration of 
Cu 2+ ions in KCUF3. Due to the splitting of the 3d states in the octahedral field within 
the CuF 6 octahedra, the hole at each magnetic Cu 2+ ion occupies one of two degenerate 
e g orbitals. The superexchange coupling (fT3~I) is usually analysed in terms of e g holes 
in this case pU], and this has become a textbook example of spin-orbital physics by 
now [TTlfTOT]. 

The bilayer spin-orbital model is obtained following [70]; it describes S = \ spins 
with the Heisenberg SU(2) interaction coupled to the e g orbital r = \ pseudospins, with 
orbital operators ([3]) obeying much lower cubic symmetry of the orbital exchange: 

(<i>ll7 L V 7 

+ (r 2 + r 4 )Hfe) Q-rW) Q - rf) } 

-E Z J2^ C K (55) 

The energy scale is given by the superexchange constant (TT5]) . with t standing here 
for the (ddo~) effective hopping element. The terms proportional to the coefficients 
{r 1 ,r 2 ,r 4 } originate from the charge excitations to the upper Hubbard band [70] which 
occur in d 9 d 9 =f± d 8 d^° processes and depend on Hund's exchange (fT6l) parameter, with 
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111 

1 — 3rj 1 — 1] 1 + rj 

Note that t- operators are not independent because they satisfy the local constraint 
= 0. The bilayer model (155]) depends thus on two parameters: (i) Hund's 
exchange coupling 77 ffTB"]) . and (it) the crystal-field splitting of e g orbitals E z j J. 

The II*. operators stand for projections of spin states on a bond (ij) on a singlet 
(Ilfj) and triplet (II^) configuration for £ = | spins, i.e., 

nf > = ( j - s, • s,) , = (| + s, • s,-) . (57) 

Their form suggests that singlet spin correlations will play an important role in particular 
parameter regimes. The usual on-site MF approximation captures only global symmetry 
breaking in the bilayer, with essentially four different magnetic phases: (z) two G-AF 
phases with FO order characterized by either x or z orbitals occupied by the holes and 
stable at large values of \E Z \, (ii) the A-AF phase stabilized by finite Hund's exchange 
rj ffT6]) near the orbital degeneracy E z — 0, and (Hi) the FM phase which has the lowest 
energy at a sufficiently large value of rj. It is clear that a better approach than the 
MF approximation with on-site order parameters (S*) and (t^ ) has to be employed 
to capture subtle effects of spin fluctuations which may stabilize VB or resonating VB 
(RVB) phases. Indeed, it has been shown that the phase diagram obtained in the on-site 
MF approach changes drastically and is much richer when the cluster MF approach is 
used instead [36], see below. 

A more sophisticated approach which goes beyond the single-site MF approximation 
takes a cubic 2x2x2 cluster as a reference, with eight corner sites coupled to 
their neighbours along the bonds in the ah planes by the MF terms. This choice is 
motivated by the form of the Hamiltonian containing different interactions in three 
different directions, and the cube is the smallest cluster which couples the ah planes and 
does not break the symmetry between the a and b axes as it contains equal numbers of 
a and b bonds. In the considered case of a bilayer there are no further neighbours of the 
cube along the c axis. After dividing the entire bilayer into identical cubes which cover 
the bilayer lattice, the Hamiltonian fl55l) can be written in a cluster MF form as follows, 

meC 

where the sum runs over the set of cubes C, with each individual cube labeled by C m G C. 
Here contains all bonds that belong to a given cube C m and the crystal-field terms 
oc E z at cube sites, i.e., it depends only on the operators on the cube sites, while 
contains all bonds outgoing from a given cube m and connecting it with neighbouring 
clusters. 

The basic idea of the cluster MF approach is to approximate "H^f by T-L^ containing 
only operators from the cube m. This can be accomplished in many different ways 
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depending on which type of symmetry breaking is investigated. A natural choice is to 
take of the following general form [36] . 



"f — a,b 
ieCm 



containing SU(2) symmetry-breaking spin field S?, orbital field and spin-orbital field 
Sfrf . Coefficients {aj, fej, cj, <fj} are the Weiss fields determined self-consistently for 
various values of the parameters {E z /J, 77}. Note that we introduce here spin-orbital 
field SjTj because, as seen a posteriori in some phases, spins and orbitals alone do not 
suffice to describe the symmetry breaking which may occur when these operators act 
together. 

The standard way, as in any MF approach, is to derive self-consistency equations 
for the Weiss fields. This can be done in a straightforward fashion: we take the operator 
products from anc ^ divide them into a part depending only on operators from the 
cube m itself, and from a neighboring cube n. This procedure can be applied to all 
operator products in ti^" and full H^f can be recovered in the MF form fl59|) . After 
repeating this procedure for all clusters, one finds a set of commuting cubes interacting 
in a self-consistent way. The Weiss fields at site i, 

«Z = ^2 + n)u] + -(r 2 - n)sj , (60) 
b l = - ( r 4 + r x )u] - ^(r 2 - n)s], (61) 
^(3n-r 4 )t7 + ^(r 2 + r 4 ), (62) 



7 



d] = - -(n + r^uju 1 ^ - ^(r 2 - r^u] + s]u^ 
-^(r 2 + r 4 )(C,-tZ) + i(r 4 -3r 1 )tZt 7 



1 1 1. 1 



Si = 


(S?) . 


t 1 ■ = 




7 

U m,i = 


(si 



- l(3n + 2r 2 + r 4 ), (63) 
are determined together with the order parameters at site i, 

(64) 
(65) 

5 -***))■ (66) 

The next crucial step is to impose a condition that {s],t],uj} are related to the 
order parameters obtained on the internal sites of the considered cluster. Thereby it is 
convenient to assume that two neighbouring cubes can differ in orbital (and spin-orbital) 
configuration by the interchange of a and b direction, i.e., 

$1 = ±Si. t] = t]n,i, uj = ±u r mi (67) 
with 7 being the complementary direction in the oh plane to 7, i.e., (7, 7) = (a, b), (b, a). 
This relation gives the same results as the one when a ■H- b symmetry in the cube is not 
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Figure 10. (Color online) The phase diagram of the cluster MF Hamiltonian 
of the cP spin-orbital model for a bilayer, with independent spin (1651) . orbital (|66[) 
and spin-orbital (|66|) MFs. VB phases (mainly PVB and VBz) with spin disorder are 
stable in the light shaded (yellow) area, and phases with spin-orbital entanglement are 
indicated by dark grey (orange) shading. This figure is reproduced from [36j . 



broken, but keeps the whole system a-b symmetric in the other case. It is also important 
that the spin-orbital field in not factorized (but surprisingly it turns out that such a 
factorization does not prevent spin-orbital entanglement to occur [36]). 

When the spin-orbital MF (|66|) is not factorized but calculated self-consistently, 
one finds the phase diagram of figure [Ull Here two G-AF phases for E z > with 
occupied x orbitals, A-AF with accompanied C-AO order and FM are familiar from the 
simplest on-site MF approach, and appear also in the 3D KK model [70]. The part of 
the diagram for 77 < 0.2 and E z j J < 0.5 is dominated by two new phases compared with 
a one-site MF approach (shown in [36]): VBz phase with z orbitals occupied stable for 
negative E z , and plaquette VB (PVB) phase in a range of E z ~ 0, both with vanishing 
magnetization. For small E z /J < 0.25 one finds that quantum fluctuations included 
within the present approach select the former spin-disordered VB2; phase. The latter 
PVB phase has singlets formed on the bonds either along the a or b direction of the 
cluster, depending on the cube. This phase breaks the a -B- b symmetry locally but 
the global symmetry is preserved thanks to the 7r/2 rotation of neighbouring clusters 
mentioned above. The orbitals take shape of cigars pointing in the direction of the 
singlets. 

The most important result presented in figure [TOl are the regions of stability of three 
new entangled phases: ESO, EPVB and PVB-AF, obtained only when the spin-orbital 
order parameter uj (1661) is not factorized into the spin and orbital part. ESO stands 
for entangled spin-orbital phase and is characterized by relatively high values of spin- 
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orbital order parameters, especially for high values of r\ when other order parameters 
are close to zero. This phase contains singlets along the bonds parallel to the c axis, 
its magnetization vanishes and orbital configuration is nonuniform. EPVB stands for 
entangled PVB phase and resembles it, but has in addition a finite spin-orbital field, and 
weak global AF order. Finally, a different type of phase with spin-orbital entanglement 
is the PVB-AF phase connecting PVB and G-AF in a smooth way, stable only if 77 is 
large enough. 

Note that with the exception of this last phase, the other two entangled phases arise 
near the quantum critical point (QCP) which is found at 77 = in the MF approach 
and moves to a finite value of 77 ~ 0.12 when quantum fluctuations on singlet bonds 
are explicitly included. Therefore, the phase diagram of figure [10] implies that singlet 
formation suppresses frustration caused by Hund's exchange coupling and moves the 
region of the most frustrated interactions to finite 77. This shows once again that the 
simple single-site MF approach is insufficient to describe faithfully the phase diagram 
of the present spin-orbital model ( )55l) . We suggest that similar entangled phases are 
expected in the 3D model which should be investigated within the cluster MF approach 
in the near future. 

5.2. Spin-orbital resonating valence-bond liquid 

In this section we consider another example of spin-orbital entangled states which are 
found in the d 1 model on the frustrated triangular lattice [371138] . as realized in the (111) 
planes of NaTiC^. In the limit of large intraorbital Coulomb interaction U intersite 
charge excitations are again transformed away and one finds the following effective 
Hamiltonian [37] , 



a bond in different processes. Here the interaction may arise either from superexchange 
via oxygen orbitals due to transitions via the effective hopping t that follows from the 
d — p hybridization, or via direct (kinetic) exchange between t2 g orbitals with flavour 
7 active on a bond along the direction 7 in a triangular lattice via the hopping t'. 
The parameter a interpolates between the superexchange (a = 0) and direct exchange 
(a = 1) limit. It is the first parameter of the present model (|68p and is given by the 
ratio of these two hopping elements: 



The superexchange involves two t2 9 orbital flavors (jH]) different from 7 on each bond that 
are not active in direct exchange. Consequently, the superexchange is more quantum, 
similar to the superexchange in titanates or vanadates, while the direct exchange is 
more classical in the orbital channel, bearing some similarity to the superexchange for 




a = 



t 2 + t' 2 ' 



(69) 
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e g orbitals analysed above in the bilayer KK model. The second parameter of the spin- 
orbital model f )68|) is Hund's exchange rj ffl6|) . as in the KK model. More details on the 
model and its derivation can be found in [37] . 

In the subsequent sections we focus first on the frustrated interactions in the d 1 
model ( Fl3l) at rj = limit, and we give here its explicit form in this case: 



Hq—J 
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(70) 



The summations include the bonds (ij) || 7 of a triangular lattice, with 7 = a,b,c 
labeling three directions. This case is rather special as the multiplet structure collapses 
to a single excitation with energy U (spin singlet and triplet excitations are then 
degenerate), and the Hamiltonian simplifies. The operators electron number 

operators for the orbital flavour 7 at site z, and nf 1 ^ is the density in the remaining two 
orbitals, which is related to by the local constraint, 



(7) 



n i7 + n\ 

The scalar products of the orbital operators in (1701) . 
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(71) 



(72) 



involve two active orbital flavours on superexchange bonds. For a bond along the axis 7 
orbital operators at site i are defined by the electron creation {a\, b\, c}} and annihilation 
{aj, 6j, Cj} operators for electrons with a given flavour. For instance, for the bonds along 
the a or b axis they are: 



T+ = b] 



T: 



ii, 
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(73) 
(74) 

(75) 



The labels 7 7^ 7 in the quantum fluctuating part oc a/«(1 — a) refer to the two orbital 
operators on each bond along the direction 7 involved in the fluctuating operators 
defined in fl73l) . Orbital fluctuations are the only processes contributing to the mixed 
exchange terms in this limit (77 = 0). 

A remarkable feature of the Hamiltonian in the limit of 77 = ( 1701 is the lack 
of higher symmetry in any of the points when a is varied. Even at a = 0.5, where 
all electron transitions have the same amplitude, no higher symmetry occurs as the 
superexchange (a = 0) and direct exchange (a = 1) result from quite distinct processes 
and involve different subsets of orbital flavours which cannot be transformed one into 




Figure 11. (Color online) Evolution of the properties of a single bond 7 = c for the 
spin-orbital model (1701) on a triangular lattice as a function of a at r\ = 0: (a) energy 
spectrum E n (solid lines) with degeneracies given by numbers; (b) spin $y (|19l) (filled 
circles), orbital TJy (fSU|) (empty circles), and spin-orbital Cy (l?T1) (x) correlations. 
The transition between the two distinct regimes occurs by a level crossing at a — 0.5. 
For a < 0.5, the two types of dimer wave function [(ss/ot) and (os/st)] are degenerate 
(d = 6) for resonating orbital configurations {ab}, while at a > 0.5, the nondegenerate 
spin singlet is supported by c orbitals occupied at both sites [(ss/cc)]. This figure is 
reproduced from [37) . 



the other. The only analytical solution was found in the a = 1 case, where at r] = 
the extremely degenerate GS is a liquid of hard-core dimers |102j . This degeneracy is 
removed at finite r] > 0, and a VB crystal with a large unit cell of 20 sites is formed. 

To understand a subtle interplay between the quantum superexchange at 77 = 
and more classical direct exchange in the limit of rj = 1 we consider first a single bond 
oriented along the c axis. In the superexchange limit the active orbitals are a and b, 
while only c orbitals contribute to the direct exchange ( ffPl . A single bond gives the GS 
energy E = — J, both in the superexchange (a = 0) and in the direct exchange (a = 1) 
limit, see figure ITTT a). These two limits differ in a fundamental way — the GS at a = 
has degeneracy d = 6 due to two complementary triply degenerate wave functions, with 
spin singlet and orbital triplet (ss/ot), and spin triplet accompanied by orbital singlet 
(st/os), while in the opposite a = 1 case limit spin singlet is accompanied by frozen FO 
order of active c orbital (ss/cc). Energy increase when the a = 0.5 point is approached 
from either side indicates increasing frustration. 
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The different character of wave functions prevents any energy gain that might result 
from orbital fluctuations in the mixed exchange term, and the GS energy E increases 
linearly toward Eq = —\J when the point a = 0.5 is approached from either side. At 
this point the interactions are maximally frustrated, degeneracy is d — 7, and a QPT 
between the two described GSs takes place. Also the remaining part of the spectrum 
(excited states) is symmetric with respect to the a = 0.5 point. Altogether, the evolution 
of the spectrum with a demonstrates not only that superexchange and direct exchange 
are physically distinct, excluding each other and unable to contribute at the same time, 
but also that the two wave functions optimal in either limit are extremely robust. 

The above interpretation is consistent with the spin, orbital, and composite spin- 
orbital correlation functions for the considered bond, defined in equations (|19p - (|2ip . In 
the entire regime of < a < 0.5, averaging over two degenerate (ss/ot) and (st/os) 
wave functions leads to equal spin and orbital correlation functions = TV,- = — ~, see 
figure [TUb). As a singlet for one quantity is matched by a triplet for the other one, the 
two sectors are strongly correlated, and one also finds — —\. Although individual 
quantum states may be written here as products of the respective spin and orbital states, 
this result suggests that strong spin-orbital entanglement in expected in the GS of a 
larger system. It may be shown that entanglement arises there mathematically because 
the GS is a resonating superposition of a number of configurations [38] which do not 
couple here with one another and form degenerate states for a single bond. Note that 
the obtained value, = — |, is the same as in the SU(4) model, see section |3T2| and 
reflects maximal possible entanglement. 

By contrast, for a > 0.5 the above degenerate states favoured by superexchange 
become excited states, and the spin-singlet GS (SV,- = — |) is the lowest, see figure fTTTb). 
The orbital configuration is characterized here by a rigid order of c orbitals, (n ic nj C ) = 1, 
which quenches all orbital fluctuations. Thus the spin and orbital parts are trivially 
decoupled, giving CV,- = 0. Finally, at the transition point a = 0.5 one finds degeneracy 
d = 7 and averaging over all degenerate states yields Sy = —0.321, TV,- = —0.214, and 
somewhat reduced composite function CV,- ~ —0.23, as one of the degenerate states 
gives here CV,- = 0. Summarizing, the regime of entangled spin-orbital states a < 0.5 is 
complemented by a factorized (disentangled) (ss/cc) wave function. 

Below we analyse how the above two distinct regimes of the d l model are modified 
in case of large number of neighbours and frustration in the triangular lattice. We 
consider two clusters with PBC: a hexagonal cluster of N = 7 sites (N7), and a rhombic 
cluster of N = 9 sites (N9). Unfortunately, the cluster with even N = 12 number of sites 
is too large and would require very extensive calculations, but the presented clusters 
are sufficient to demonstrate essential differences between the nature of interactions in 
both limits. Due to the PBC all sites and bonds are equivalent, so bond correlations 
are all the same and independent of the bond direction 7. Each t<i g orbital is occupied 
on average by 1/3 electron, but similar to a single bond the GSs are manifestly different 
in the limits of a = and a = 1, see the insets in figure fT2Tb). 
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Figure 12. Bond correlations for clusters with PBC at r/ = — spin S (|19|l. orbital 
T (|20p. and spin-orbital C f|2 1 1) : (a) hexagonal cluster of N = 7 sites, and (b) rhombic 
cluster of N = 9 sites. Vertical lines indicate the exactly determined range of possible 
values of each quantity at a given a that follows from the GS degeneracy. The insets in 
(b) indicate typical orbital patterns in the superexchange (a — 0) and direct exchange 
(a = 1) limit for the rhombic N9 cluster; a more symmetric hexagonal cluster N7 
is obtained by removing top-right and bottom-left sites. These data were presented 
in (38]. 

In case of N7 cluster all three directions are equivalent and the spin correlations are 
AF and constant, Sij ~ —0.11 independently of a. Note that this value can be deduced 
as follows: In the low-spin phase one has St — § total spin, and one can determine the 
intersite spin correlations using the following identity [38J: 

S 2 = 7S 2 i + A2(S i -S j ). (76) 

Due to PBC every pair of sites stands for a nearest neighbour bond and this implies 
the above relation. The value (S{ ■ Sj) = —3/28 = —0.107 obtained from it reflects 
AF spin correlations but is much reduced from the classical limit of —0.25 found in 
the Neel state on a square lattice, in agreement with high frustration of the triangular 
lattice. Note that the present situation is radically different from the ID SU(4) model, 
where frustration is absent but spin correlations are reduced by their coupling to the 
orbital correlations. In the N9 cluster spin correlations are in a similar range but vary 
somewhat with a as the cluster shape breaks the symmetry between bonds in three 
nonequivalent directions {a,b,c}. These correlations are weaker (<S ~ —0.090) at a = 
than in the N7 cluster and become more pronounced (S ~ —0.144) when a ~ 0.6, while 
orbital fluctuations gradually weaken to C ~ —0.050 at a = 1. 
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In both clusters orbital correlations are negative, T < 0, and the Goodenough- 
Kanamori rule stating that these correlations should be complementary to spin ones 
is violated. The orbital correlations weaken in both clusters when a increases towards 
a = 1 and the superexchange interactions are less important, particularly in the N9 
cluster. Joint spin-orbital correlations are also similar in both clusters (e.g. C ~ —0.070 
at a = 0) and \C\ gradually decreases when spin and orbitals disentangle approaching 
a — 1. 

An important question is whether spin order and excitations could be described 
by an effective spin model derived from the spin-orbital model (fl~3|) . Below we show 
that this is not the case in the d 1 model on the triangular lattice. MF procedure used 
frequently leads here to [38] : 

"Hmf = ^2 {(^J 7) ) ^ ' ^ j ~ (^ 7) ) ' ^)} 

<«>I|7 

+ E{4 7) (^-4)+4 7) }- (77) 

(y}||7 

In this way spin and orbital degrees of freedom are disentangled and the model reduces to 
a superposition of the spin model with self-consistently determined orbital correlations, 
and the orbital model, with self-consistently derived spin correlations, similar to the 
decoupling of spin and orbital degrees of freedom introduced for the ID spin-orbital 
chain in section 14.31 Following the spin model, one obtains the MF spin interactions 
for N7 and N9 clusters by averaging the orbital operator (its explicit form is given 
in [3H]) over the MF GS |$o) (03) • Note that the orbital fluctuations in the term 
oc a/«(1 — a) in (!70|) contribute here as well as they couple different components of 
|$o)- 111 contrast, the exact exchange constant J exa ct ( ED} is found when the exact GS 
|$) obtained after Lanczos diagonalization is used for a given cluster. 

In figure [13] we compare the phase diagrams obtained from the above MF procedure 
and from Lanczos diagonalization for the N7 cluster. Consider first a QPT from the 
low-spin (S t = 1/2) disordered phase to the high-spin (S t = 7/2) FM phase which 
occurs for sufficiently large r\. When spin and orbital operators are disentangled ( 1771) . 
see figure [TBT a). and C = 0, it coincides with the sign change of the MF exchange 
constant Jmf ( ]T4l) and no intermediate phase (with 1/2 < S t < 7/2) is found, as in a 
spin system. 

Comparing the values of Jmf CHI) and J exa ct (|30]) found from the MF and from 
exact diagonalization of the N7 cluster, see figure [131 one finds that J exa ct > Jmf 
in a broad range of a. Therefore, the MF approximation turns out to be rather 
unrealistic and overestimates (underestimates) the stability of states with FM (AF) 
spin correlations. The value of Jmf decreases with increasing rj, but positive values 
Jmf > are found at r\ = only if 0.27 < a < 1. This demonstrates that FM 
states: (i) are favoured when joint spin-orbital fluctuations are suppressed, and (ii) are 
stabilized by orbital fluctuations close to a = even in absence of Hund's exchange. 
The transition from the low-spin (St = 1/2) to the high-spin (St = 7/2) state occurs 
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Figure 13. (Color online) Contour plots of the effective exchange constant Jmf as 
obtained for the hexagonal cluster of N = 7 sites with PBC: (a) within the MF 
calculation which includes orbital fluctuations (TT4|) . and (b) using the exact GS found in 
exact diagonalization pop . In case (a) the transition from low-spin to high-spin phase 
occurs when the exchange constant Jmf changes sign and becomes negative. Thick 
lines in (b) indicate the phase boundaries between phases with increasing intermediate 
total spin value St =3/2, 5/2 between <S t = 1/2 and 7/2 for increasing rj. This figure 
is reproduced from [38 . 

in Lanczos diagonalization at a much higher finite value of rj rj 0.14, with only weak 
dependence on a, see figure ITBTb). In addition, one finds two phases with intermediate 
spin values St = 3/2,5/2 in a range of rj values near a = 0.8. Note that the exchange 
constant J ex act changes discontinuously at the onset of the FM phase. 

We have found that the qualitative trends presented here for the N7 cluster are 
similar to those observed for the N7 cluster |103] and thus they may be considered 
generic for the present d 1 model on the triangular lattice. In both cases one finds 
that: (i) the FM phase is stable in the MF approximation close to a = and becomes 
degenerate with the low-spin phase at a = 1, (ii) the MF procedure is exact in the 
regime of FM phase, and (to) the transition to the FM phase occurs gradually through 
intermediate values of total spin St (except at a = 1). This suggests that partially 
polarized FM phase should occur in the thermodynamic limit. It arises due to spin- 
orbital entanglement which is gradually suppressed when r] increases. 
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We argue that the results presented in [3S] and [103 J provide evidence that the 
present d 1 spin-orbital model realizes a paradigm of a spin-orbital liquid phase in the 
superexchange regime, and the order-out-of-disorder mechanism does not occur when 
the Hilbert space contains coupled spin and orbital sectors. Previous search for this 
quantum state of matter in other systems, particularly in LiNi02 where e g orbitals 
are active on the d 7 configurations of Ni 3+ ions on the triangular lattice |104j . were 
unsuccessful |105[I106] . After considering the present model in more detail we suggest 
that the triple degeneracy of ti g orbitals plays a crucial role in the onset of a spin-orbital 
liquid, as the number of orbital flavours fits to the geometry of the triangular lattice. 
In contrast, the direct exchange regime is dominated by VB states with spin singlets 
accompanied by static configurations of directional orbitals providing the energy gain for 
the direct exchange. Note that the frustrated triangular lattice plays here an important 
role and removes any kind of orbital order in the entire range of a. 

6. Hole propagation in a Mott insulator with coupled spin-orbital order 

Finally, we give an example of a single hole doped into the half-filled Mott insulator with 
a pseudo-entangled AF/AO order [34J. This type of order is realized in the ab planes of 
LaV0 3 with S = 1 spins [21] and in Sr 2 V0 4 with S = ^ spins |107j . It might appear as 
entangled as it contradicts the Goodenough-Kanamori rules and has certain similarity 
to the entangled spin-orbital states in a ID chain, see section 13.31 It is well known that 
a hole doped into a quantum antiferromagnet couples to the collective (delocalized) spin 
excitations (magnons), and propagates through the lattice surrounded by a "cloud" of 
magnons which is the essence of a QP behaviour [108] • Thereby the energy scale of the 
" coherent" hole propagation is strongly renormalized from the hopping t and is given by 
the AF superexchange constant J. This QP is frequently called a spin polaron [28J and 
was observed in the photoemission spectra of the parent compounds of high-T c cuprates, 
as e.g. in Sr 2 Cu0 2 Cl 2 [T09] . 

A hole doped in a Mott insulator with orbital order could behave in a similar 
way when orbital fluctuations or interorbital hopping [30] would also repair the defects 
created by a hole. Orbital excitations are decoupled from spin ones in disentangled 
states, such as FM spin and AO order in the ab planes of LaMn03 [12]. It has been 
shown [30] that a hole introduced into such a state indeed does not disturb the FM spin 
order and couples to the collective excitations of the AO state (orbitons). Here again a 
QP is formed which is called this time an orbital polaron [110] . Usually, however, the 
bandwidth of an orbital polaron is much smaller than that of the spin one. In fact, the 
orbitons are in general considerably less mobile than the magnons (or even immobile) 
due to almost directional Ising-like superexchange [i~2H3] . Actually, one can understand 
the hole motion in this case in terms of the string picture |111[I112] . 

The orbital order created by ti g orbitals in a 2D square lattice is more robust than 
that of e g ones while only ID hopping for each orbital flavour is allowed in the former 
case [31]. The case of a hole doped into the plane with FM spin order accompanied 
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Figure 14. Artist's view of a single hole introduced into the spin and orbitally ordered 
ab plane of LaVC>3 with coexisting AF/AO order. One electron at each V 3+ ion 
occupies the xy orbital (not shown) and the other one either yz or zx degenerate 
orbital, forming the classical AO state (the projections of these orbitals onto the ab 
plane are shown) whereas their spins S = g alternate on two sublattices, forming the 
AF Neel state. This figure is reproduced from [34) . 



by the t 2g AO order, which could correspond not only to the hole introduced into the 
ordered ground state of Sr 2 V04 with t 2g orbitals but also, surprisingly, to K 2 CuF 4 or 
Cs 2 Ag0 4 with active e g orbitals. As the GS has AO order, also here a QP might formed 
due to the dressing of a hole by the collective orbital excitations. However, due to the 
specific t 2g orbital symmetries the orbitons are not mobile at all, and the QP acquires 
a finite bandwidth only due to the frequently neglected three-site terms [31 J . Thus, 
the string picture determines the nature of the t 2g orbital polarons even more than in 
systems with e g orbital degrees of freedom. 

A more complex situation arises, however, when a doped hole may couple in a t 2g 
system both to magnon and to orbiton excitations [66] • This situation is found in the 
ah plane of LaVOs, shown schematically in figure [14] and studied in [34]. The coexisting 
AF/AO order consists of the AF order induced by the superexchange interactions 
between the occupied c orbitals, and the AO order of active {a, b} orbitals. The GS in an 
undoped system is therefore a perfect two-sublattice order, with quantum fluctuations 
in the spin channel but a classical AO order in the orbital channel as the superexchange 
interations in this part are Ising-like. The coexistence of the AO and AF order is 
extremely rare as it formally violates [T] the Goodenough-Kanamori rules predicting 
complementary spin and orbital order in the GS. This GS has long-range order, with an 
up (down) spin component accompanied by a zx {yz) orbital, but this composite nature 
of the GS has no dramatic consequences as long as the system is undoped. 

The problem of a single hole doped into the ab plane of LaV03, see figure [HI 
is challenging as the hole doping occurs in the orbital {a, b} doublet and disturbs 
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locally both AF and AO order [113] . Unlike in the FM/AO case [22], for the present 
coexisting AF/AO order neither the spin nor the orbital background is transparent for 
a propagating hole, and the hole has to couple simultaneously to both magnons and 
orbitons when it moves by one lattice spacing. It is then unclear whether a QP may 
form and which role is played individually by spins and orbitals in possible formation 
of a spin-orbital polaron. 

The 2D spin-orbital t-J model may be seen as a generalization of the spin t-J 
model [61] and the orbital t2 9 t-J model [31] to the spin-orbital superexchange — it 
consists of three terms 



H = H t + Hj + H 3s , (78) 

where the last one, Hs s , stands for the three-site effective hopping. The second term 
Hj is the spin-orbital superexchange model for the -RVO3 perovskites (Tl3"]) introduced 
above in section 14.11 The first term in (178]) describes the hopping of {a, b} electrons in 
the constrained Hilbert space, i.e., in the space with singly occupied (at hole position) 
or doubly occupied sites. This means that electrons in the ah plane, which is under 
consideration here, can hop only along the b (a) direction when they carry a (b) orbital 
flavour. The c orbitals do not participate in hopping processes as they are always 
occupied by one electron ( 129]) . see [80 11114] . Hence, we arrive at the kinetic oc t part of 
Hamiltonian (1781) . 

H ^ = - t Y, V + «U i+ b CT + H.c.) V. (79) 

i,cr 

Here the constrained operators a[ a , mean that the hopping is allowed only in the 
restricted Hilbert space with not more than one {a, b} electron at each site i (a stands 
for the spin component opposite to a). Besides, since Hund's exchange coupling is large 
(Jh t) [115] . we project the final states resulting from the electron hopping onto the 
high-spin states, which occurs due to the V operators in ( 1791) . More details can be found 
in [32]. 

Low energy excitations are magnons and orbitons, with their energies determined 
by the respective exchange constants derived from the spin-orbital superexchange [3] 



_ 1-377 - V , m 

Js - 4(1 -377X1 + 2,7) J ' (80) 

Jo " (1-377X1 + 277) J - (81) 
One finds that Jo > and Js > in the expected range of 77 < 0.2, which means that 
the classical GS has indeed coexisting AF and AO order. By rotating first spins and 
orbitals at sublattice A to the FM/FO order and introducing Schwinger bosons f\ a for 
spin and t\ for orbital excitations, one finds the energies of magnons and orbitons using 
the LSW theory. As usually in a quantum antiferromagnet, magnons are dispersive and 
have a Goldstone mode, 



u± = J 8 zSJl-il, (82) 
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where 7k is given by (j53p . S = 1 and z = 4 is the coordination number. In contrast, 
the orbital excitations at energy Jq are local and have no dispersion as the orbital 
interactions are Ising-like. 

A crucial step in deriving the hole spectral function A(k, u) is the following 
representation of the electron operators in terms of the {fi a ,t\ a } (a = a,b) Schwinger 
bosons: 

4 = , (83) 

The above equations demonstrate that a spin excitation is always generated together 
with an orbital excitation, which implies that the diagrams contributing to the self- 
energy contain only vertices with two outgoing or incoming excitation lines. The factor 
follows from spin algebra for spins S = 1 in the spin-orbital model [H3J- Here the 
projection onto the high-spin states has already been done, so the projection operators 
V in H t are no longer needed. 

The spectral function can be obtained from the Green's functions which are defined 
separately for A and B sublattice and depend on the corresponding self-energy S a (k, uj) 
(a = a, b): 

G a (k,u)= 1 (85) 

to + e a (k) - E a (k, oj) 

Here e a (k) = 2rcos(2fc y ) and £ft(k) = 2rcos(2/c x ), and r = jj describes the three- 
site hopping responsible for the ID dispersion within each sublattice. In what follows 
we show only the result for r = for more clarity. Green's functions are solved 
self-consistently together with the self-energies derived in the self-consistent Born 
approximation (SCBA). After solving the Green's functions, the spectral functions for 
a hole created in a = {a, b} orbital are: 

A a (k,u) = --lim lmG a (k,u + i5) . (86) 

At t = the spectral function does not depend on the orbital flavour, and we define 
A(k,u) = A y (k,u). 

The spectral functions A(k, u) ( !86l) shown in figure [15] are representative for the 
situation in the strongly correlated transition metal oxides. They are obtained by solving 
the SCBA equations for J = OAt on a mesh of 16 x 16 k-points. Similar data for other 
values of J = 0.2t and 0.6t are presented in [3lj. A realistic value of Hund's exchange in 
LaV0 3 is r] = 0.15 [71]. The peak in the low-energy part of the spectrum, see figure [ToTc). 
has no dispersion which indicates hole confinement (a rather weak ID dispersion along 
the k x (k y ) direction for holes doped into the b (a) orbitals is obtained at finite r = \ J). 
The spin-orbital spectral functions form ladder-like spectra and have many similarities 
with the spectra obtained for the t<i g orbital t-J model [21] • However, there are few 
subtle differences with the orbital model which demonstrate that the spin-orbital case 
is more complex and the spin part also contributes. 
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Figure 15. Spectral functions for a single hole in the AF/AO phase along the T — M 
direction of the Brillouin zone for: (a) toy orbital model (|88p (top); (b) toy spin model 
dSTJ (middle); (c) the full spin-orbital model (|T3|) (bottom). Parameters: J = OAt and 
rj = 0.15 (i.e., J s = OMt and J = 0.16t). Broadening 5 = O.Olt in the definition 
of the spectral function is assumed. Image courtesy of Krzysztof Wohlfeld. 

To understand better the nature of the obtained spectra we present also the spectral 
functions obtained for the related spin and orbital problem in figure IToTa-b). These two 
models read as follows [34J: 

h s = - tj2 ptfLty + H - c -)^ + j s E s i • s j » ( 87 ) 

Ho = - 1 J2 (^l^i+a + al«i+b + H.c.) + Jo T ? T i- ( 88 ) 

Here spin operators {Si} stand for S = 1 spins, Tf are zth components of pseudospin 
T = 1/2, and the operators V project onto the high-spin states. The constrained 
electron operators, c\ a = c\ a (l — n^) in ( IH71) and b\ = b\(l — ni a ) and a\ = a\(l — nib) 
in (IHHl) exclude double occupancies from the Hilbert space in each case, similar to ( J79|) . 
The superexchange energy scale is Js (IHUl) for the spin model (IH7|) and Jo ( IHTj) for the 
orbital one ( 1881) . which mimics the influence of the orbital part on the spins and of the 
spin part (with the AF order) on the orbitals. With the present parameters one finds 
the spin-only exchange constant Js = 0.06t (somewhat higher than that deduced from 
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the observed value of the Neel temperature T N ~ 143 K in LaV0 3 [21]), and Jo = 0.16t. 

The spectral functions for the above models were obtained using the SCBA on a 
mesh of 16 x 16 k-points, following the derivations presented in [28] and [32] in the 
case of the spin and orbital model. One only has to make the following substitutions 
in the respective SCBA equations: J — > — Js, S — > 1 and (due to the quantum double 
exchange factor) also t — > tj y/2 in the spin case [28], and J — > — Jo, E — > 0, and r — > 
in the orbital case [32J. The spectral functions obtained for the above toy models are 
quite different from those obtained for the spin-orbital t-J model. The spectral function 
for the spin model, see figure [ToT b). consists of a QP and the incoherent part at higher 
energies, both having considerable k-dependence. It is evident that this k-dependence is 
removed by the orbital interactions. The orbital model at r = has a strictly localized 
spectral function with a ladder spectrum [28], see figure IToTa). The spin-orbital spectral 
function, see figure fToTc). resembles qualitatively the ladder spectrum found for the 
orbital model in figure [ToTa). The momentum dependence of both the QP state and the 
incoherent part of the spectrum obtained for the spin model, see figure [15T b). is entirely 
suppressed. These results demonstrate that the quantum spin fluctuations are to a large 
extent quenched in the spin-orbital model by the simultaneous coupling of the hole to 
both spin and orbital excitations. It is remarkable that the hole still couples to the spin 
degrees of freedom by generating string potential due to defects created by hole motion. 
Thus, the string which acts on the hole moving in the plane with AF/AO order is of 
the composite spin-orbital character. This not only explains the peculiar correspondence 
between the orbital and spin-orbital model, but also shows that the spins play an active 
role in the lightly doped spin-orbital system. 

We remark that the orbitally induced string formation considered here could be 
understood as a topological effect. It happens even if the orbital excitation energy is 
turned to zero, i.e., when the hole moves in the orbital sector incoherently. Hence, the 
mere presence of orbitals is sufficient to obtain the (almost) classical behaviour of a hole 
doped into the GS with AF/AO order. This result suggests that further investigation of 
the hole propagation in spin-orbital systems is a fascinating subject for future studies. 

7. Discussion and summary 

Spin-orbital entanglement discussed in this topical review concerns entanglement on 
the exchange bonds, similar to entanglement in e.g. spin singlets which build VB 
states [39lH0]. We have shown that this concept is important is several spin-orbital 
models and in general either the ground state or excited states are entangled. Two 
recent examples of entangled ground states were presented: the Kugel-Khomskii model 
on a bilayer, and the spin-orbital d l model on the triangular lattice, in addition to the 
well known ID SU(2)®SU(2) model. It may be expected that more spin-orbital model 
systems with entangled ground states will be found in the near future. Whether or not 
such states are indeed realized in nature depends on the coupling to the lattice. For 
instance, the e g orbital order is robust in LaMnOs [12] and KCuF 3 [99], and for this 
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reason the ground states of these compounds are disentangled. 

We have also demonstrated that the Goodenough-Kanamori rules are violated in 
the regime of weak Hund's exchange in several situations. It is in this regime that the 
mean field decoupling procedure of spin and orbital operators fails and the magnetic 
properties can be determined only by solving the full entangled spin-orbital many-body 
problem. Also for weak (or vanishing) Hund's exchange spin- orbital liquid phase is 
stabilized by spin-orbital entanglement in the d l spin-orbital model on the triangular 
lattice. 

However, a frustrated lattice does not guarantee that a disordered spin-orbital 
liquid-like state arises. A good counterexample is the frustrated lattice of alkali RO2 
hyperoxides (with R=K, Rb, Cs), where the interactions induced by the lattices compete 
with the superexchange and stabilize the orbital order, lifting the geometric frustration 
of the lattice [116] . Then the spin-orbital interactions may be consideed as disentangled, 
and the Goodenough-Kanamori rules require certain reinterpretation. In fact, they have 
been generalized by including large interorbital hopping terms |116j . 

It has been shown that ground states which contradict the celebrated Goodenough- 
Kanamori rules may appear in two different situations: (i) either when quantum 
fluctuations in the spin-orbital system are strong and stabilize the entangled ground 
state with coexisting AF and AO order, or (ii) when the AF superexchange interactions 
follow from charge excitations in other orbitals not involved the actual orbital order, 
or (Hi) when the coupling to the lattice stabilizes the AO order in the regime where 
the AF spin correlations are expected due to superexchange. In these latter situation 
the ground state is more classical and not entangled as it follows from interactions 
induced by the lattice and not from the spin-orbital superexchange. Also in the regime 
of large Hund's exchange Jh, ground states with coexisting FM and FO order may 
appear |117j . similar to the FM/FO states in the SU(2)xSU(2) model. They contradict 
again the Goodenough-Kanamori rules, but are in fact disentangled — here quantum 
fluctuations are absent and play no role for their stability. 

On the example of the -RVO3 perovskites we have shown that the experimental 
data in this unique family of correlated oxides indicate that low-energy excited states 
are entangled and the energies of spin and orbital excitations are similar. This happens 
because lattice distortions and interactions of t2 g orbitals with Jahn- Teller modes are 
weak and only the electronic interactions such as superexchange decide about the system 
behaviour. In the -RVO3 systems the temperature dependence of the optical spectral 
weight, the phase diagram as a function of the ionic radius of R ions, and the dimerized 
magnon excitations, are all determined by the presence of entangled states with low 
excitation energies. 

The interplay between spin, orbital and spin-orbital excitations poses a very 
interesting problem for future theoretical studies as well as an experimental challenge. 
We also point out that in case of FO order accompanied by AF interactions, spin 
fluctuations couple to orbital excitations and thus they cannot be considered separately 
|118j . in spite of formal separation of spin and orbital degrees of freedom in the 
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ground state. Another important finding is recent observation that composite spin- 
orbital excitations fractionalize in the quasi-lD Mott insulator Sr 2 Cu03 |119j . The 
nature of spin-orbital excitations and the circumstances of their possible decay are very 
challenging and unresolved questions. The case considered here of a hole moving in 
the AF/AO order |Mj seems to suggest that, at least in doped systems, spin-orbital 
excitations may play a very important role in transport as they may not decay and 
impose some topological constraints on carrier propagation. 

The present topical review focused on spin-orbital entanglement on the 
superexchange bonds, but a different kind of entanglement arises in presence of on- 
site relativistic spin-orbit interaction. In this case a Kramers doublet gives the lowest 
energy states of a single ion and these states are next considered to derive interactions 
between neighbouring ions |120j . Such local on-site entangled states play an important 
role and decide about the magnetic properties of Sr2ir04 with weak ferromagnetism and 
Sr 2 V04 with hidden spin-orbital order. In a ID superexchange model with Ising-like 
orbital superexchange locally entangled states introduce orbital dynamics and lead to a 
phase diagram with a novel phase having long-range Neel order of spin and orbital 
angular momenta [121] . Excitations is such models are a challenge and they are 
under investigation at present. If localized states with an effective angular momentum 
Jj = 2Si — Li at site i (with S and L being spin and orbital operators) are considered, 
one finds a gapped spectrum at finite spin-orbit coupling [122J. 

Finally, we would like to emphasize that a better understanding of the concept of 
spin-orbital entanglement is important and could help to make progress in other fields as 
this subject is interdisciplinary. Recently a scheme to generate spin-orbit-path hybrid 
Greenberg-Horne-Zeilinger entanglement [123] was proposed [124] for photons which 
are entangled with different degrees of freedom. In quantum chemistry entanglement is 
considered in chemical bonds which are classified using measures of electron correlation 
[125] and entanglement [126] . Spin-orbital-like entanglement is also applicable to nuclear 
systems, where nucleons possess as well two degrees of freedom — spin and isospin. 
It has been established that the entanglement length of the nucleons is significantly 
larger than that one expects [78J. Entanglement may also play an important role in 
quantum computations if the spin state used for information storage would be measured 
by investigating orbital qubits in entangled states [127] . 
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